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Abstract 



We study the Gross-Pitaevskii (GP) energy functional for a fast rotating Bose-Einstein condensate 
on the unit disc in two dimensions. Writing the couphng parameter as 1/e^ we consider the asymptotic 
regime e — J- with the angular velocity fi proportional to (e^| log£|)~^. We prove that if fi = 
no(e^| loge|)~^ and ilo > 2(37r)~^ then a minimizer of the GP energy functional has no zeros in an 
'•f/^ annulus at the boundary of the disc that contains the bulk of the mass. The vorticity resides in a 

complementary 'hole' around the center where the density is vanishingly small. Moreover, we prove a 
lower bound to the ground state energy that matches, up to small errors, the upper bound obtained 
. . from an optimal giant vortex trial function, and also that the winding number of a GP minimizer 

^ around the disc is in accord with the phase of this trial function. 

-J—^ 
^ MSC: 35Q55,47J30,76M23. PACS: 03.75.Hh, 47.32.-y, 47.37.+q. 
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1 Introduction 

An ultracold Bose gas in a magneto-optical trap exhibits remarkable phenomena when the trap is set 
in rotational motion. In the ground state the gas is a superfluid and responds to the rotation by the 
creation of quantized vortices whose number increases with the angular velocity. The literature on this 
subject is vast but there exist some excellent reviews [XI ICoi iFtel] . The mathematical analysis is usually 
carried out in the framework of the (time independent) Gross-Pitacvskii (GP) equation that has been 
derived from the many-body quantum mechanical Hamiltonian in [LSYj for the non-rotating case and in 
[LSj for a rotating system at fixed coupling constant and rotational velocity. When these parameters tend 
to infinity the leading order approximation was established in [BCPYj but it is still an open problem to 
derive the full GP equation rigorously in this regime. The present paper, however, is not concerned with 
the many-body problem and the GP description will be assumed throughout. 

When studying gases in rapid rotation a distinction has to be made between the case of harmonic traps, 
where the confining external potential increases quadratically with the distance from the rotation axis, 
and anharmonic traps where the confinement is stronger. In the former case there is a limiting angular 
velocity beyond which the centrifugal forces drive the gas out of the trap. In the latter case the angular 
velocity can in principle be as large as one pleases |Fe2| . Using variational arguments it was predicted in 
[FBj that at sufficiently large angular velocities in such traps a phase transition changing the character of 
the wave function takes place: Vortices disappear from the bulk of the density and the vorticity resides 
in a 'giant vortex' in the central region of the trap where the density is very low. This phenomenon was 
also discussed and studied numerically in both earlier and later works jFJSi IFZ[ IKBl IKFl IKTU| , but 
proving theorems that firmly establish it mathematically has remained quite challenging. 

The emergence of a giant vortex state at large angular velocity with the interaction parameter fixed 
was proved in |R1| for a quadratic plus quartic trap potential. The present paper is concerned with 
the combined effects of a large angular velocity and a large interaction parameter on the distribution 
of vorticity in an anharmonic trap. In particular, we shall establish rigorous estimates on the relation 
between the interaction strength and the angular velocity required for creating a giant vortex. The 
physical regime of interest and hence the mathematical problem treated is rather different from that in 
[Rlj but a common feature is that the annular shape of the condensate is created by the centrifugal 
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forces at fast rotation. By contrast, the paper |AAB| focuses on a situation where, even at slow rotation 
speeds, the condensate has a fixed annular shape due to the choice of the trapping potential. The basic 
methodology of the present paper is related to that of |AABj and jIMll IIM2J but with important novel 
aspects that will become apparent in the sequel. 

As in the papers |CDY11 ICY| . that deal mainly with rotational velocities below the giant vortex 
transition, the mathematical model we consider is that of a two-dimensional 'flat' trap. The angular 
velocity vector points in the direction orthogonal to the plane and the energy functional in the rotating 
frame of reference i^ 

(1.1) 



{I 



f'-^i*] := / df<^ |V^r-2r2-^*L* + £"^ 



*r} 



where we have denoted the physical angular velocity by 251, the integral domain is the unit disc B = 
{r £ K^ : r < 1}, and L = —if A V is the angular momentum operator. It is also useful to write the 
functional in the form 



£ 



GP 






df 



V-iA * 



-n' 



m 



s-^W^ 



(1.2) 



where the vector potential A is given by 



A-~Q Af^ fire 



i5- 



(1.3) 



Here (r, "(9) are two-dimensional polar coordinates and e^ a unit vector in the angular direction. The 
complex valued \1/ is normalized in L^{B) and the ground state energy is defined as 



E' 



GP 



inf s^'^m. 

11*11, =1 



(1.4) 



The minimization in (1.4) leads to Neumann boundary conditions on dB. Alternatively we could have 



imposed Dirichlet boundary conditions, or considered the case of a homogeneous trapping potential as 
in |CDY2j . Our general strategy applies to these cases too, but with nontrivial modifications and new 
aspects that are dealt with in separate papers |CPRY] . 
We denote by vJ/GP a minimizer of (1.2), i.e., any normalized function such that f'^^[^'^^] = E'^^ . 



The minimizer is in general not unique because vortices can break the rotational symmetry |CDY1[ 
Proposition 2.2] but any minimizer satisfies the variational equation (GP equation) 



- AvI/GP - 2d • L vI/GP + 2e-2 |vi/GP|^ xpGP ^ ^ 
with Neumann boundary conditions and the chemical potential 

^GP.^^GP^l /dflvI/GP'^ 



GP^GP 



(1.5) 



(1.6) 



The subsequent analysis concerns the asymptotic behavior of vJ/GP ^nd E'^p as £ — ?> with ft tending to 
oo in a definite way. 



As discussed in |CDY1] the centrifugal term — 51^r^|^p in (1.2) creates for 51 > e ^ a 'hole' around 



the center where the density |v[/G!P|2 jg gxponentially small while the mass is concentrated in an annulus 
of thickness ~ e51 at the boundary. Moreover, by establishing upper and lower bounds on E'^^ , it was 
shown in |CY) that in the asymptotic parameter regime 



log£|«51«(e2|log£|)- 



(1.7) 



^The notation "a := b" means that a is by definition equal to b. 
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the ground state energy is to subleading order correctly reproduced by the energy of a trial function 
exhibiting a lattice of vortices reaching all the way to the boundary of the disc. It can even be shown that 
in the whole regime (1.7 1 the vorticity of a true minimizer ^'-'^ in the annulus is uniformly distributccrl 
with density Q/n. 

1.1 Main Results 

In this paper we investigate the case 

" = ^^ (1-8) 

e-^l loge| 

with J7o a fixed constant and prove that for Qq sufficiently large a phase transition takes place: Vortices 
disappear from the bulk of the density and all vorticity is contained in a hole where the density is low. 

To state this result precisely, we first recall from [CDYlj that the leading order in the asymptotic 
expansion of the GP ground state energy is given by the minimization of the 'Thomas- Fermi' (TF) 
functional obtained from £'^'^ by simply neglecting the kinetic term: 

£TF[^] .^ f ^^ {^n\^p + e-^p^} , (1.9) 

Jb 

where p plays the role of the density |^p. The minimizing normalized density, denoted by p"^^, can be 
computed explicitly (see Appendix). For ft > 2{y/7Te)^^ it vanishes for r < i?h with 

l-Rl = {2/^){en)-^ ^ n^^e\ loge], (1.10) 

while for R\^ < ?' < 1 it is given by 

p^^{r) = \e^n\r'-Rl). (1.11) 

Our proof of the disappearance of vortices is based on estimates that require the TF density to be 
sufficiently large. For this reason we consider an annulus 

A:^{r : Ri, + e\\oge\~^ <r<l} (1.12) 

where p"^^{r) > 17o£~^| loge|~'^. Our result on the disappearance of vortices from the essential support 
of the density is as follows: 

Theorem 1.1 (Absence of vortices in the bulk). 



If Vt is given by (1.8) with flo > 2{3tt) ^, then '^'^^ does not have any zero in A for e > small enough. 



More precisely, for f £ A, 

\\^^^{r)\' - pTF(,)| < c^ « p^^ir). (1.13) 

Remark 1.1 (Bulk of the condensate). 

The annulus A contains the bulk of the mass. Indeed, because fj p^^ = 1 — o(l) and ^>^^ is normalized, 

( |1.13[ ) implies that also /_^ |*c;P|2 ^i_ ^^^^y 



The proof of Theorem 1.1 is based on precise energy estimates that involve a comparison of the energy 
of the restriction of '3>*^p to A with the energy of a giant vortex trial function. By the latter wc mean a 
function of the form 

$(7=) = /(?^exp(ir2i9) (1.14) 



•^ Cf. |CYI Theorem 3.3]. This theorem is stated only for the case O < e ^ but it holds in fact true in the whole region 
(TtI i [R3] . 
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with / real valued and f2 G N. It is convenient to writaj 



(1.15) 



with w € Z and use w as label for the J7-dependent quantities in the sequel. We define a functional of / 

by 

SS^'lf] := £°P[/exp(zJ7i?)] = / df {|V/p + {n' r'^ - 2nCl),f + e-^f} . (1.16) 



Minimizing the functional £^^[f] (see Proposition 2.3), that is convex in /^, over L^ -normalized / for 
fixed ft, we obtain an energy denoted by E^^ and a minimizer that is rotationally symmetric, without a 
zero for r > 0, and unique up to sigrjj We denote the unique positive minimizer hy g^. If < oj < Ce~^ 
the corresponding density g^ is close to p^^ (see Remark 1.3 below). 

It is clear that 

SGP < £GP[^^ exp(^^!^?)] = ^^^ (1.17) 

for any value of uj and hence 

^GP < ^GP .^ -^f ^GP^ (1^18) 



Our second main result is a lower bound that matches (1.181 up to small errors (see Remark 1.4 
below) and an estimate of the phase that optimizes E^^ . 

Theorem 1.2 (Ground state energy). 

For rtf) > 2(37r)^^ and e > small enough the ground state energy is 



E 



GP 



E 



GP 



O 



I log e 1 3/2 



£l/2(log|log£|)2 



Moreover E'^^ = E!^^ , with Wopt € N satisfying 

2 



w, 



opt 



3v7re 



(l + 0(|loge| 



-1/2- 



(1.19) 



(1.20) 



Because ^'^^ does not have any zeros in the annulus A the winding number (degree) of ^p'^^ around 
the unit disc is well defined. Our third main result is that the giant vortex phase [fl] — ojopt is, up to 
possible small errors, equal to this winding number. 



Theorem 1.3 (Degree of a minimizer). 

For J7o > 2(37r)^^ and e > small enough, the winding number of '$ is 

deg {*°P, as} = [O] - c^opt(l + o(l)) 



(1.21) 



with u! opt o,s in (1.201 



The following remarks are intended to elucidate the phase difference Wopt and to justify the claim that 



the remainder in (1.19) is, indeed, a small correction to E'^'^ . 



Remark 1.2 (Giant vortex phase). 

For a fixed trial function / the energy (1.16) is minimal for 



n = n[ dfr-'f +0(1) 



(1.22) 



^We use the notation ' [ • ] ' for the integer part of a real number. 

^Instead of requiring / to be real we could have required / to be radial in which case the minimizer is unique up to a 
constant phase. 
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If /^ is close to p^^ and thus essentially concentrated in an annulus of width 
follows that 

o< {n-Ci)^o{e-^). 



(en) = 0(e|log£|), it 
(1.23) 



The leading term 2/(3-y/7r)£ ^ in (1.20) can, indeed, be computed from (1.221 with /^ 

simply, for /^(r) increasing linearly from at r = i?h to its maximal value at r = 1. It is worth noting 



TF 

p , or more 



that if /^(r) were constant in the interval [i?h,l], the leading term for uj according to (1.221 would be 
2/{-s/tt£) and the phase Ct therefore approximately equal to (il/Tr) x area of the hole. The optimal Cl is 
closer to U, because of the inhomogeneity of the density in the annulus. 

Remark 1.3 [Giant vortex density functional). 



The functional (1.16) can also be written 



Jb 



(1.24) 



with 



B^{r) -.^nr-n/r. 



(1.25) 

li oj = 0{e^^), then B^{r) = 0{e^^) close to the boundary of the disc. Neglecting the term B^f^ as well 
as the gradient term in ( 1.24[ ) in comparison with the centrifugal and interaction terms leads to the TF 
functional evaluated at f''. This makes plausible the assertion above that g^ is close to p'^^ if cj = 0{e^^) 
and this will indeed be proved in Section [2] 

Remark I.4 (Composition of the ground state energy). 



If one drops the gradient term in (1.24) but retains the term with B^, one obtains a modified (and 
w-dependent) TF functional 



Jb 



P- 



-2„2 



pI 



(1.26) 



IB 

Its minimizer and minimizing energy E'^^ can be computed explicitly (see the Appendix) and one sees 
that the energy is minimal for oj = 2/(3v'7r£)(l + o(l)). Denoting the minimal value by E"^^ , we have 

^TF ^ ^TF ^ 0(^-2) 



(1.27) 



while 



E^^ = -n^ - 4/{3VTr)e-^n = 0{e-^\ logej-^) + 0{e~^\ loge 



(1.28) 



The term 0{e~^) in (1.27) is the angular kinetic energy corresponding to the third term in (1.24). The 
difference between E*^^ and E^^ is the radial kinetic energy of the order (efi)^! logej — 0{e^^\oge\~^). 
The remainder in (1.19) is thus much smaller than all terms in E'^^ . It is also smaller than the energy a 
vortex in A would have, which is | loge| x density — ©(e^^). 



1.2 Proof Strategy 

We now explain the general strategy for the proof of the main results. 

The first step is an energy splitting as in LM\ We write a generic, normalized wave function as 
^(r) = gu:{'''')w{r) with a complex valued function w. The identity 



/df|V*|2= /dr(-Ag^)g^|u;|2+ j Argl\Vw\ 

JB JB JB 



(1.29) 
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holds by partial integration and the boundary condition on g^j . Using the variational equation for g^^ and 

(1.30) 



writing w{'r) — exp(if2'i?)z;(r) and fi as in (1.151, we obtain a splitting of the energy functional 



with 



where 



£u[v] 



2\2 



dr gi\Vv\' - 2gtB^ ■ [iv^v) + ^(1 - \v\') 



BUr) ■■= B^{r)e^ = pr - {[n] - uj) r"!] 



and we have used the notation 

{iv,Wv) := ^i{vWv* — w*Vw). 

Since the divergence of the two-dimensional vector field g^B^^ vanishes we can write 

2glB^ - V^^„ 
with a scalar potential F^ and the dual gradient V^ = {—dy,dx)- Stokes theorem then gives 



£M 



df{gl\Vv\^+F^cuT\iiv,\7v)}~ / daF^{iv,drv)+ / df^{l-\v\^) 



OB 



|2\2 



(1.31) 

(1.32) 
(1.33) 

(1.34) 
(1.35) 



Here and in the rest of the paper 'curl (if, Vw)' stands for the 3-component of V A (iw, Vw), d-r is the 
tangential derivative and dtr the Lebesgue measure on the circle, i.e, given a ball Bh of radius R centered 
at the origin, dr '■= R~^d§ and dcr := Rd'd on dBR. 

The giant vortex trial function 5^ exp(if2'(9) gives an upper bound to the energy and hence we see 



from (1.30) that 



£M<^ 



(1.36) 



where u is the remaining factor of vl/^^ afte r the giant vortex trial function has been extracted, e.g., 
\1/*^^ = g^j exp(if2i9) u. The proof of Theorem 1.1 is based on a lower bound on £uj[u] that, for fig large 



enough, would be positive and hence contradict ( 1.36 ) if u had zeros in A. The main steps are as follows: 



• Concentration of the density: There is an annulus A^ slightly larger than A^ so that the supremum 
of l-i-CPp over B\A ij^ 0(g°° ). The same holds for gl provided w < Ce'^ This implies that for 
V = u the integrations in (|1.31 ) can be restricted to A up to errors that are 0{e°°). 



• Optimization of the phase: Choice oiui = 0{e ^) so that F^^, chosen to vanish on the inner boundary 
of ^, is small on dB. 

• Boundary estimate: Using the variational equations for vj/GP q^^i^ g^ QT^id the sniallness of F^ on 



the boundary it is proved that the boundary term in (1.35) is small. This step uses the Neumann 



boundary conditions on dB in a strong way and its extension to the case of Dirichlet boundary 
conditions is a nontrivial open problerq^ 

• Vortex halls: Isolation of the possible zeros of u in 'vortex balls' and an estimate of the first term 
in ( 1.35 1 from below in terms of the infimum of g^ on the vortex balls and the winding numbers of 



u around the centers. For the construction of vortex balls a suitable upper bound on the last term 



in (1.35) is essential. 



^We use the symbol 0{e°°) to denote a remainder which is smaller than any positive power of e, e.g., exponentially 
small. 

®This problem has recently been solved |CPRY) . 
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• Jacobian estimate: Approximation of the integral of F^^ against the vorticity curl(jM, Vw) by a 
sum of the values of Fi^ at the centers of the vortex balls multiplied by the corresponding winding 
numbers. 



• Gradient estimate: Estimates on Vu leading to a lower bound on the last term in (1.351 that, for 
Qq large enough, excludes zeros of u. 

A key point to notice is that the vortex ball construction in combination with the jacobian estimate leads 
to a cost function defined as 

HUr):^yi{r)\loge\+FUr). (1.37) 

The first term is, to leading approximation, the kinetic energy of a vortex of unit strength at radius r, 
the second term is the gain due to the potential energy of the vortex in the field 2g^Bi^. If H^^ is negative 
at some point, the energy may be lowered by inserting a vortex at this point. A positive value of H,^ 
means that the cost of the kinetic energy outweighs the gain. Note that i/„ depends on VIq through F^; 
in fact it is not difficult to see thatjl] 



\F.{T)\<^^9l{.r)\\oge\. 



(1.38) 



Hence, if VLq is large enough, the cost function is positive everywhere on the annulus and vortices are 
energetically unfavorable. An upper bound, 2(37r)~^, on the critical value for Q,q is computed in the 
Appendix. This upper bound is in fact optimal, as demonstrated in [R2] . The proof requires additional 
ingredients. 

The construction of vortex balls is a technique introduced independently by Jerrard . Jli and Sandier 
[Saj in Ginzburg-Landau (GL) theory, whereas the jacobian estimate originates from the work |JS| . Both 
techniques are described in details in the monograph [SS]. This method has been applied in |AAB) and 
|IM1[ |IM2] to functionals that at first sight look exactly like (1.31). There is, however, an essential 
difference: While in |AABj and [IMll IIM2| the size of the relevant integral domain is fixed, the weight 
function g^ is in our case concentrated in an annulus whose width tends to zero as e — >■ 0. This shrinking 
of the width of the annulus and the large gradient of g^ are the main reasons for the complications 
encountered in the proof of our main theorems. 

To appreciate the difficulty the following consideration is helpful. As already noted, the concentration 
of g^ on the annulus A of width 0{e\ loge|) implies that B^ 



0{e 



on this annulus. From (1.361 and 



the Cauchy-Schwarz inequality applied to the second term in (1.31), one obtains the bound 



dr%(l 



\u\'?< 



C 



(1.39) 



Now, although g'^ ^ (£^)^ ^ e "^l logej ^ on A, the estimate (1.39) is not sufficient for the construction 



of vortex balls on A. In fact, ( 1.39 ) is compatible with the vanishing of u on a ball of radius 0{e\ loge|), 
i.e., comparable to the width of A, while for a construction of vortex balls one must be able to isolate 
the possible zeros of u in balls of much smaller radius. 

The solution of this problem, elaborated in Sectio n [4} involves a division of the annulus into cells of 
area ~ (e^)^ with the upshot that a local version of (1.39) holds in every cell. The local version, 



/ ■ 

J Cell 



df^-§{l~\u\r< 



"^Y < ^ y- (number of of ceUs) ^ 



(1.40) 



means that, in the cell, u can only vanish in a region of area 0{£^\ logep) that is much smaller than the 
area of the cell, i.e., ^ e^| logep. Hence the vortex ball technique applies in the cells where (1.40), or a 



^Here, as in the rest of the paper, C stands for a finite, positive constant whose value may vary from one formula to 
another. 
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sufficiently close approximation to it, holds. The gist of the proof of a global lower bound, that eventually 
leads to Theorems |1.1| and |1.2[ is a stepwise increase in the number of cells where an estimate close to 



(1.40) is valid until all potential zeros in the annulus are included in vortex balls. The proof of Theorem 
|1.3| involves in addition an estimate on the winding number of u. 

In the sketch of the proof strategy above we have for simplicity deviated slightly from the actual 
procedure that will be followed in the sequel. Namely, for technical reasons, we find it necessary to 
restrict the considerations to a problem on the annulus before the splitting of the GP energy functional 
as in (1.301. This means that, instead of the function guj^^t ^^d the optimal phase Wopt defined above, we 



shall work with corresponding quantities for a functional like (1.24 1 but with the integration restricted 
to A. The main reason for this complication is lack of precise information about the behavior of \I/^p 
in a neighborhood of the origin. In fact, the distribution of the 'giant vorticity' of ^"^^ in the hole is 
unknown. In particular it is not known whether ^*^^ vanishes at the origin like g^j, which has a zero 
there of order Cl. 

1.3 Heuristic Considerations 

A heuristic argument for the giant vortex transition, based on an analogy with an electrostatic problem, 
has been given in |CYj and goes as follows: Writing a wave function as ^{f) ^ |^(r)| exp{i(p{f)) with a 



real phase ip the kinetic energy term in (1.2) is 



IdrllVl^lp + lM/piv^p-Ip} (1.41) 



with A{r) = ilre^. We focus on the second term and consider the situation where the phase (p contains, 
besides the giant vortex, also a possible contribution from a vortex of unit degree at a point vq in the 
annulus of thickness ~ (^£)~ • The phase can thus be written 

ip{r) := (I'd + arg(f - fo). (1.42) 

where a,Tg{x,y) = arctan(2//a;). The modulus |^| vanishes at tq and is small in a disc ('vortex core') of 
radius ~ \/ s/Vl ^ £'^"| loge|^" around ro. The phase arg(r — rg) can be regarded as the imaginary part 
of the complex logarithm if we identify M? with C The corresponding real part, i.e., conjugate harmonic 
function, is log \r — ro| which is the two-dimensional electrostatic potential of a unit point charge localized 
in tq. Likewise, the conjugate harmonic function of Cl-d is filogr, i.e., the electrostatic potential of a 
charge U, placed at the origin. By the Cauchy-Riemann equations for the complex logarithm we can write 

V(^ = dr^p Br + r^^d^ip e^ = drX &i9 ^ '''^^d^X eV (1-43) 

with 

X(r) :=f^logr + log|f-fo|. (1.44) 

After a rotation by 7r/2, i.e., replacement of e^ by eV at every point, the vector potential also has an 
electrostatic interpretation: Qrer is the electrostatic field of a uniform charge distribution with charge 



density fi/Tr. Employing (1.43) we now have 

\Vip-A\^ = \Vx-^rer\'^ (1.45) 

and E(f) :~ Vx — ^JreV is the electric field generated by the point charges and the uniform background. 
We can now apply Newton's theorem to argue that the effect of the giant vortex, i.e., the first term in 



(1.44), is to neutralize in the annulus the field generated by the uniform charge distribution in the 'hole'. 



i.e., the second term in ( 1.45 1, provided fi is chosen to match the area of the hole times the uniform charge 
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density. By the same argument the point charge at tq neutrahzes, outside a disc of radius ^ 1/v ri, the 
effect of one unit of the continuous charge. 

The total charge of the continuous distribution in the annulus is Q ^ fl{ile)^^ ^ e^^. Inserting the 
vortex reduces the effective charge by one unit so the corresponding energy gain is ^ e^^. On the other 
hand, the energy associated with the electrostatic field from the point charge outside the vortex core 
(that is cut off by the modulus of the wave function) is ~ |^P| loge| ~ (^e)| loge|. The condition for the 
cost outweighing the gain is ile\ logej > e~^, i.e., 



1 

e^ I log i 



^^T^ITZZZ^^ (1-46) 



marking the transition to the giant vortex phase. 

It is also instructive to consider a version of the functional ( 1.31 ) where the variables have been scaled 
so that the width of the annulus becomes 0(1). We define 

i:={en)-'^n^'e\\oge\, s:=r^r, g{s) ■= 9u>m, v{s) := viis), B{s) := W^iis). (1.47) 



Then (1.31) is equal to 



£[v] := I dsf ||Vt)p - 2B • (w, Vv) + ^(1 - 1*1')') (1-48) 

with B :~ Bg-i a ball of radius £^^ ~ J7o s^^\ logej^^. For the sake of a heuristic consideration we now 
assume that g^ is constant, ^ £^^ , on the annulus of width £ and zero otherwise. We define a new small 
parameter 

e -.^ e£-'/^ ^ e'/^\loge\-'/^ (1.49) 

and note that 

\B\ ^ e-H ^ n-W logel + 0(log | loge|)). (1.50) 

We are thus led to consider the functional 

ds||Vup-2B-(u',Vi)) + 4(l- \v?Y\ (1-51) 

on an annulus A of width 0(1) and an effective vector potential of strength ©(rig^^l logej). This is 
reminiscent of the situation considered in |AAB| and jIMli IIM2] where the domain is fixed and the 
rotational velocity is proportional to the logarithm of the small parameter. Moreover, increasing f2o 
decreases the coefficient in front of the logarithm. Hence for large J7o this coefficient is small and if 
the analysis of jAABj and |IMH IIM2J would apply, one could conclude that vortices are absent. This 
reduction of the problem to known results is, however, too simplistic because the annulus A is not fixed: 
Although its width stays constant, its diameter and hence the area increases as e~^|loge|~^. A new 
ingredient is needed, and in our approach this is the division of the annulus into cells as mentioned in the 



previous subsection. In the scaled version (1.481 of the energy functional these cells are (essentially) of 
fixed size. When writing the actual proofs we prefer to use the original unsealed variables but the picture 
provided by the scaling is still helpful, in particular for comparison with |AAB| and |IMHITM2] . 

Remark 1.5 (Alternative approach). 

After the submission of this paper we learned [J3| of a possible alternative approach to prove the absence 
of vortices in the bulk, relying on |j2l Lemma 8] (see also |AJR[ Lemma 4.1]). This method could replace 
some of the arguments in Section |4] and lead to a shorter proof of our Theorem |1.1| but is likely to yield 
worse remainder terms in the energy (Theorem 1.2). We also stress that the tools developed in Section 111 



of the present paper are an essential input in the paper iR2) , where it is proved that vortices do appear 
in the bulk if fio < 2(37r)^^. For this Lemma 8 in [J2] would not be sufficient. 
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1.4 Organization of the Paper 

The paper is organized as follows. In Section [2] we gather some definitions and notation that are to 
be used in the rest of the paper. We also prove useful estimates on the matter densities, both for the 
actual GP minimizer [vj/CPp ^^^ ^^^ ^^le 'giant vortex profiles' g^. In Section^we introduce the auxiliary 
problem that we are going to study on the annulus A and prove that it indeed captures the main energetic 
features of the full problem. Section |4] is devoted to the study of the auxiliary problem via tools from the 
Ginzburg-Landau theory. We conclude the proof of our main results in Section Is] Finally an Appendix 
gathers important facts about the TF functionals that we use in our analysis, as well as the analysis of 
the cost function that leads to an upper bound on the critical speed. 

2 Useful Estimates 

In this section we state some useful estimates which are going to be used in the rest of the paper. Some of 
them are simple consequences of energy considerations and in particular energy upper bounds, whereas 
others depend in a crucial way on the variational equations solved by ^'^^, etc., and apply therefore only 
to (global) minimizers. 

In the first part of the section we investigate the properties of any GP minimizer: Starting from 
simple energy estimates, we prove the concentration of \1/^^ on an annulus of width 0{e\ loge|) close to 
the boundary of the trap and in particular its exponential smallness inside the hole. This is the major 
result about ^^^. 

The second part of the section is devoted to the analysis of the densities associated with the giant vortex 
energy. We shall introduce first some notation and then, exploiting energy considerations, discuss the 
properties of those minimizers: Most of them are very similar to the ones proven for ^ , as, e.g., the 
exponential smallness, but we also need other general properties as, for instance, an a priori bound on 
the gradient of the densities and an estimate of their difference from the TF density p^^ . 

2.1 Estimates for GP Minimizers 



We briefly recall the main notation: Given the GP energy functional £'^^[^] defined in (1.2), we denote 
by E'^^ its infimum over L^-normalized wave functions. Any GP minimizer is denoted by vp^^ and solves 



the variational equation (1.51 



The TF functional i£"'"^[p] was introduced in (1.9) and E'^^ and p^^ stand for its ground state energy 

and density respectively (see also the Appendix). 

Any further label to the functionals £'~'^ and £^^ , as, e.g., £§^ , denotes a restriction of the integration in 

the functional to the domain V. The same convention is used for the corresponding ground state energies 

and minimizers. 

Finally we use the notation B{f, g) for a ball of radius g centered at f, whereas Br is a ball with radius 

R centered at the origin. 

The starting point is a simple GP energy upper bound proven for instance in [CY;, Proposition 4.2], 
i.e., 

A straightforward consequence of such an upper bound is that the L^ norm of vJ/GP jg concentrated in 
the support A"^^ := {r : i?h < r < 1} of the TF minimizer. At the same time the bound implies a useful 
upper bound on j'!/ |: 

Proposition 2.1 (Preliminary estimates for ^I**^^). 
As e ->-0, 

l|l*^"P-P^'IL^(B) <0(1)' ll*^'ll'==(B) < VIl^.b, (i + 0(v/^1^)) . (2.2) 
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Proof. In order to prove the first statement it is sufficient to use the fact that p^^ is the positive part of 
the function (e^/2)(/i'^^ + fl^r^) together with the normaUzation of p^^ and ^^^ and estimate 



dr [|vI/GP|2_pTF] 



■'"''rA + VWl-^ I df |V|/GP|%TF < 



■^ 



GPI 



„TF| 



* lU + ll^^'lb-^ / ^r\^>- \- p' 

JB 

] - eV^^ - e'^' / dfr^l^pGPp ^ ^2 (^tf [|^gP|2] _ ^tf^ < 

£2 (i^GP _ ^TF) < ^^ (2.3) 



by (2.1) 



The proof of the second inequahty is similar to the proof of Lemma 5.1 in |CY| and involves the 
variational equation (1.5). We define 

pGP ,^ |,l;GP|2^ (2.4) 

The crucial point is that at any maximum point of p^^ in the closed ball B, it has to be Ap*^^ < 0. This 
is trivially true in the open ball B but can be extended to the boundary thanks to Neumann boundary 
conditions: Since 9^^*^^ = at the boundary 96, which implies drp'^^ = there, p*^^ can have a 
maximum at Fq S dB only if V/9'^^(ro) — and therefore Ap'^^(ro) < 0. 
From the variational equation (1.5) solved by ^'^^, one can estimate 



ApGP < 4e-2 (e^^GP ^ ^2^2^2 _ 2pGP^ „gp 



-2 /'^2,,GP I ^2r)2 



GP\ „GP 



[e~p"^ + s~\rr- --zp"^ ) p"^ S 4e ~[e~p,"^ +eVr~-2p^^)p 
by using the properties 

-ApGP ^ „^GP*^^GP „ ^GP^^GP* _ 2 |V*«P|' , 

Now since Ap^^ < at any maximum point of p*^^ , one immediately has 



(2.5) 



(2.6) 



,GP| 



-2|,,GP 



On the other hand 



„TF II /^„ri /^„— 1| 1 — , „l — 1 



,TF| 



P'\l) 



Cen^Ce-'\\ose\ 



(2.7) 



and the difference between the chemical potentials (see (1.6) and (A.2) for the definitions) can be esti- 
mated as follows 



,GP ,,TF| 



< E 



GP 



TF 



E 



Ce~ 



TF 



1 1/2 



+ * 



GP| 



TF||l/2 



|,I,GP|2_^TF|| < 



GP ,,TF|\l/2 



l+e3|loge||^GP_^TF|) 
Ce-5/2| logepi/^ (l + e^/^l loge|i/2|^iGP - pTF|i/2 



which yields 



and thus the result. 



|^GP_^TF|<^^-5/2|i^g^|-l/2 



(2.8) 

n 



A consequence of the L^ estimate in (2.2 1 is that the L^ norm of v&'^p is concentrated inside the 
support of p^^ or, in other words, the mass of ^'^^ inside the hole Br^^ is small. Indeed one can easily 
realize that the first inequality in (2.2 1 implies 



* 



GP| 



\L^BaJ 



<0{^/^\\^\), 



(2.9) 
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since, due to the normalization of p and its support in A 



iTF 



^ 



GPI 



\L^BnJ 



* 



GPI 



IL2(^TF) 



* 



GPI 



\LHBnJ 



dr 



A"^ 



(|vl/GP| 



TF 



1 > 



l*''''llL^(6«J + l-^(v/^^Og^)' 



where in the last step the Cauchy-Schwarz inequality and the size 0{e\ loge|) of A"^^ is used. 

This simple estimate can be refined (see, e.g., |CDY11 Proposition 2.5]) and one can actually show 
that ^'^^ is exponentially small in e inside Br^^. The next Proposition is devoted to the proof of two 
pointwise inequalities of this kind: 

Proposition 2.2 (Exponential smallness of '3>*^p inside the hole). 

As e — > and for any fCzB, 



*°P(f)| < Ce-^\\oge\-^ exp 



1-r^ 



(2.10) 



Moreover there exists a strictly positive constant c such that for any 0(e^' ^) < r < Rh — 0(e^' ^), 

|*^^(r)|'<C£-i|logeri exp{-^}. (2.11) 



Remark 2.1 (Comparison between (2.10) and (2.11)j. 

At first sight the pointwise estimates in (2.101 and ( |2.11[ ) look very similar. In fact, since 1 — i?^^ = 
0{e\ loge|), one can easily realize that the first one yields a much better upper bound than the latter as 
soon as 1 — r^ ^ ©(e^/^), i.e., in particular for l — r^= 0(1). The main drawback of the first inequality is 
however that it becomes much weaker and even useless if one gets closer to the radius it!h- More precisely 
as soon as 1 — r^ < 0{£) the bound is no longer exponentially small in e. On the opposite the second 
inequality has no r dependence and a worse coefficient in the exponential function but it holds true and 
yields some exponential smallness up to a distance of order e^^^ from the boundary i9Sflj^ . Because of 
this fact the second inequality will be crucial in the reduction of the original problem to another one on 
an annulus close to the support of the TF minimizer. The first inequality, on the contrary, will be useful 
in a region far from the boundary of the hole ^BIi^^ and close to the origin, where the second bound 
provides a worse estimate. 



Note also that the factor in front of the exponential in both (2.10) and (2.11) is essentially given by 



the sup of |\|/'^Pp over the domain B. In the case of the first inequality this is actually needed in order 
to make it meaningful on the whole of B. 



Proof. The starting point of the proof of (2.101 is the inequality (2.5) together with the estimate (2.8), 
which yield 

(2.12) 



^TF 



Ap^^<4e-^ p^'^ (r) + C£-5 1 ioge|- 



^GP 



P 



GP 



where we have set 



pTF(,) 



[A^^^ + f^V], 



(2.13) 



which coincides with the TF density p^^ inside A^^ and is negative everywhere else. More precisely, for 
any r such that r^ < R'^ — e, one has 



Ap' 



GP 



2e-^|loge|~V^^ <0, 



(2.14) 
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since in that region 



-TF 



p''{r) + Ce-^\\oge\-^ < -£^0^(1 - o(l)) <- 



2£|log£|2' 



On the other hand the function 



W{r) := exp 



r^ -1 



satisfies for any r < 1 

-AW + 2e~^\loge\'^W > 



C 



£2 I log£P 



-£|log£| 



-e-^] W>0. 



If we then multiply W by ||p'^^||oo we get then a supersolution to (2.14), so that by the maximum principle 
(see, e.g., [EJ Theorem 1 at p. 508]) 



* (?1 < * 



GPI 



Wir), 



for any r < Rh — e. However W{r) is an increasing function and W{^yR'^ — e) = ||p'^^||oo by construction, 
which implies that the estimate trivially holds true for any r Cz B. The estimate (2.101 is a consequence 



of (2.2 1 and the inequahty 1 — i?^ = £| log£| :^ e. 



Concerning the proof of the second inequality (2.11 ), we first notice that (2.12 1 implies that, for any 
'e Bu,_^ such that r < Ry, - 0{e'^/'^), 



since in that region 



^TF 



p-(r) = 



ApGP<-C£-"/«|log£|-V''P, 



{r' Rl) < 



(2.15) 



(2.16) 



£5/6|log£|2' 

For any tq such that rp < Rh — C'(£^^^) we can thus consider the annulus X = {r £ B : r £ [ro — 9,i"q + q]} 
for some g ~ 0(£^/^), such that r^ + g < R\^ — C'(£^/^), i.e., the outer boundary of the annulus is still at 
a distance of order £^/^ from dBfi^^ , and it is straightforward to verify that the function 



U{r) 



„GP| 



exp 



(r-ro)^ 



r5/2 



(2.17) 



satisfies 

-AC/ 



4(r-ro)2 



2 ^ 2(r-ro) 



U{r) > -C 



r5/2 



U{r) > 



Ce-^/'^Uir) > -C£-i^/S| log£|-'C/(r), (2.18) 



where we have used the fact that |r — ro| < g inside I and g = 0{e'^^^). Denoting now V{'r) = 
pGP(^jfj _ U{r), we aim at proving that V^ < 0, but one has, for any r € X, 

-AV{r) < -C£-i^/^|log£|-V(r). 

Now at any maximum point of V in the interior of I it must be AV < 0, which implies V < because of the 
above inequality. Hence it remains only to prove that y < at boundary dl since the function might have 
a positive maximum there. However by construction y < at dX since U{ro — g) = U{ro + g) = ||p'^^||oo 
and thus V{r) is negative everywhere. In particular p (r) < U{r) for all f such that r — vq and 



0(£'^/6), which yields (|2.11|. 



D 
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2.2 The Giant Vortex Densities 



We now discuss some basic properties of the energy functional ( 1.16 ) together with those of an analogous 
functional where the integration is restricted to an annulus A = {r : i?< < r < 1} C ^B with _R< < i?h 
and |i?< — i?h| ^ £| loge]- A precise choice for i?< will be made at the beginning of Section IS] (see (3.2 1 ) 
but the results contained in this section apply to any i?< satisfying the above conditions. 



As indicated at the end of Section 1.2 the restricted functional is for technical reasons actually more 



useful for the proof of the main results than the original functional. It is defined in analogy with (1.161 



and (1.24 1 for real valued functions / on A by 



cGP 



f] 



A 



df {|V/p + m - Lofr-^f - 2nm - uj)f + e-2/4} 



^ dr {|V/|^ - n\'f + Blf + s-'f} , (2.19) 



with B^{r) given by (1.25). Wc use y^ as a subindex to label quantities associated with the annulus and 
thus denote the infimum of the functional (2.19) by 



pGP 



inf 

11/112 = 1:/ = / 



cGP 

^A,uj 



[/] 



(2.20) 



while E^ stands for the corresponding infimum of (1.16). 



We shall use the short-hand notation f^^, E'^^ ^ g^, etc., for quantities related either to (1.16) or (2.19) 
e.g., a statement about Qi, is meant to apply to both g^^ and gA.uj- 

Proposition 2.3 (Minimization of £^^). 

For any uj ^TL such that |aj| < 0{e~^), the ground state energies E^ satisfy the estimates 



^TF 



^GP 



^TF 



,TF 



E'' <E^' <E'f'' <E^' +Oie-'\loge\-')<E'' +Oie-') 



(2.21) 



The minimizers g^ and gA.u exist, are radially symmetric and unique up to a global sign, which can be 
chosen so that they are given by positive functions solving the variational equation 



Ag* 



m^uf 



-2„3 



2n{[n] -uj)g^ + 2e^^g'l 



^ GP 

M* 9*, 



(2.22) 



where the chemical potentials ji2^ are fixed by the L^ normalization of gi,, i.e., ft^^ — E^^ + e^^||(7j,|||. 
In addition gi, are smooth and increasing and satisfy Neumann conditions at the boundary OB, i.e., 
9^3^(1) = 0; gA,uj satisfies an identical condition at the inner boundary as well, i.e., drgA.uj{R<) = 0. 

Proof. The lower bounds to the ground state energies are simply obtained by neglecting positive terms 
(the kinetic energies) in the functionals: In the case of E^^ (the other case is identical), one has 



pTF 



Sl^''[f]>E''[f']>E''>E 

where we refer to the Appendix for the simple proof of the last inequality. 

The upper bound can be easily obtained by testing the functionals on suitable regularizations of 



(see (A.6) for its definition and (CY| . where such regularizations are performed on y p^^): The main 
correction to the energy is due to the radial kinetic energy of the regularization and one can easily realize 
that this energy can be made of order 0(£~^| loge|~^) times | loge|. Indeed, \/~p^ is a monotone function 
going from to ©(e"^/^] loge|~^/^) in an interval of size ©(ejlogel) and if it were smooth its kinetic 
energy would hence be ©(e"^] loge|~^). The exctra factor |log£| is due to the regularization close to 



Giant Vortex - CRY- December 6th, 2010 



16 



The last inequality in (2.21 ) is also discussed in the Appendix but it is basically due to the estimate 

1-8^(^)1 <n\r-'^ -r\+C\uj\+0{l) < Ce-\ (2.23) 

for any a; such that |a;| < 0{e^^) and f ^ A. 

We remark that the restriction of the integration to the annulus A has no effect because the support 

of the trial function can be assumed to be contained inside the support of p^^ , i.e., the region where 



r > R^. On the other hand by (A.9), R^ > i?h - 0(e^| loge|^), for any w e Z such that |a;| < Ce" 



which implies that the support of p^^ is contained inside A. 

Existence and uniqueness of the minimizers trivially follow from strict convexity of the functionals 
with respect to the density /^, which can be made clearer by writing them as 



4^^[Vp]=/dr{|VV^r 



+ {[n] - ufr-^p + e-2^2| _ 2r!([r!] - u), 



(2.24) 



where p = f'^ and we have used the L^ normalization of /. Similarly the radial symmetry of 5* can be 
proven by averaging over the angular variable and exploiting the convexity of the functional. 



All the other properties of g^,, including the variational equations (2.22), are trivial consequences of the 



minimization: Positivity can be proven by noticing that the minimizers g^ are actually ground states 
of suitable one-dimensional Schrodinger operators and therefore cannot vanish except at the origin (see, 
e.g., |LL1 Theorem 11.8]). Smoothness follows from (2.22) by a simple bootstrap argument, etc. 
The only property which requires a brief discussion is the monotonicity and we state it in a separate 
Lemma. D 



Lemma 2.1 (Monotonicity of the density). 

Let p{r) > be the L^ -normalized minimizer of 



dr r 



d^P 



dr 



+ ^P + bp^ 



in H^{B \ Br) C\ L'^(B\ Br) with a, 5 > and < R < 1. Then p is monotonously increasing in r. 

Proof. After a transformation of variables r'^ n- s and considering p as a function of s, the functional 
takes the form 



ds 



R^ 



dyp 



ds 



+ -p + bp' 



and the normalization condition is 



ds p = const. 



fl2 



Suppose the assertion is false. Then p has a maximum at s = si for some R < si < 1 and a local 
minimum at some S2 with si < S2 < 1. For < e < p(si) — p(s2) consider the set Ig = {s < S2 : 
p(si) — e < p{s) < p{si)}. Then, because p is continuous. 



$(e) 



ds p{s) 



is strictly positive and $(e) — > as e — )■ 0. Likewise, for (5 > we consider Js — {s > si : p{s2) < p{s) < 
p{s2) + S} and the function 

T{5):= [ dsp(s) 
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that has the same properties as $. 

Since both $ and F are continuous, strictly positive and tend to zero when their arguments tends to zero, 

there exist e, ^ > such that p{s2) + 5 < p{si) — e and <i>(e) — T{S): Because p is continuous one can, for 

any given ^ > 0, always find an e > such that the equality <i>{e) = r(^) holds and £ — > as (5 -> 0. The 

inequality p{s2) +S< p{si) — e is fulfilled if S, and hence also e, are small enough. 

We now define a new function p by putting 



Pis) :- 



We now compute the energy of the new density p. Note that it belongs by definition to the minimization 
domain and it is normalized in L^ because $(e) = r(5). 

The kinetic energy of p vanishes in the intervals Tg and J7j and the potential term for p is strictly smaller 
than for p because 1/s is strictly decreasing and the value of p on Xg is larger than on J7j. Finally 

1 fi 

2 ^ A„ „2 



P{si) - e. 


if s e Zj, 


P{S2)+S, 


if s e Js, 


p{s), 


otherwise 



ds p^ < / ds p 

because when modifying p to define p, mass is moved from Ig to J^s where the density is lower. 
Altogether the functional evaluated on p is strictly smaller than on p and this contradicts the assumption 
that p is a minimizer. D 

As next we compare the densities gl with the TF density and prove exponential smallness in the hole. 
The analogue of Proposition |2.1| is 



Proposition 2.4 (Preliminary estimates for g*). 
As e — > and for any lu £ "E such that \lu\ < 0{e^^), 



„TF| 



\lhb) 



Oil), 



2 



(B) 



< 






(B) 



(l + 0(7^1^). (2.25) 



Proof. The proof of the L estimate is exactly the same as for the first inequality in (2.2). The only 



E ) < C because of (2.21). Here the condition on w, i.e.. 



U! 



difference occurs in the energy remainder on the r.h.s. of (2.3) which can now be bounded by e^(£'p^ 



< Ce ^, is used as in Proposition 2.3 



The sup estimate can be proven by applying the same argument used to prove the second inequality 

^, one obtains the inequality 



in (2.2) to the variational equations (2.22) solved by 5*, exploiting as well Neumann boundary conditions: 

uj)'^r 



Indeed thanks to the monotonicity of the potential ([fJ] 
-Agl<Ae-'ie'i,^^ + 2s'niM- 



■)-e^iM-Lo 
which as in the proof of Proposition |2.1| implies 

e^p^^ +2e^ni[n]-Lu)-e^i[n]-Lj)^ 



29I 



9l 



\9i. 



< 



(2.26) 



The difference between the chemical potential can be estimated as in (2.8) 
sothat||5.||L<P^^(l)(l + o(l))- 



(2.27) 

n 



As for v[/G;p i-jj^g above estimates imply the exponential smallness of the densities g^, inside the hole. 



Moreover the L^ estimate (2.2) can be refined to a pointwise estimate inside A^^ and one can actually 



prove that there is a region where the difference between g-^ and p"^^ is much smaller than p"^^: 
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Proposition 2.5 (Exponential smallness of 5* inside the hole). 

As e — > and for any f (z B 



9l{r)<Ce-^\loger^exp 



l-r^ 



Moreover there exists a strictly positive constant c such that for any r < i?h — ©(e^'^), 

gl{r)<Ce-'\\oge\-'exp[--^}. 



Proof. The proof can be easily reduced to the proof of Proposition 2.2 by using (2.271. 



Proposition 2.6 (Pointwise estimate for g*). 
As £ —i' and for any w G Z with \uj\ < 0{e^^) 

laUr) ~ P^"W| < Cs'\loge\'ir' - Kr'" p"^ {r) 

for any f G A^^ such that r > i?h + 0(e^/^| logep). 



(2.28) 
(2.29) 

n 

(2.30) 



Proof. The proof is similar to the proof of Proposition 1 in |AAB] but for the sake of completeness we 
bring here all the details. 



The variational equation (2.22 1 can be rewritten in the following form 



^9* = -^ [p*{r)- gl] g*, 



where the function pi, is given by 



1 /^2,-,GP I ^2o2^2 ^2 n2 



p^{r):^l{s'f,r+e'n 



Blir)) 



Moreover by (2.231 and (2.27) 



„TF| 



|£,oc(^TF) 



< e^ lA^P - M^n + '^(l) < 0{e-^/^\ logeri/2). 



(2.31) 



(2.32) 



(2.33) 



On the other hand, for any r <E A^^ such that r — Rh > 0{e^/^\ logep), p^^{r) > 0{e ^^^) and therefore 

Mr) > P^^(0(1 - C\ logej-i/^) > Ce-i/2 (2.34) 

and in particular p^ is strictly positive for such r, which is crucial in order to apply the maximum principle. 



The pointwise estimates are indeed proven by providing local super- and subsolutions to (2.31 ) 



For the upper bound, we consider an interval [ro — S,ro + 5], where _Rh + C'£'^/^| logep + 5 < ro < 1 — S 
with ^ ^ 1, and the function 



W{r) := v^p^ro + 5) coth 



coth ^ 



p4i) 



\r -ro\ 



P*{ra + 6) 



3fc 



v/2p4ro+<5) 



One has (see |AS| Proof of Proposition 2.1]) for any r e [ro — 6, tq + S] 



AVF > ^ {p^ro + 6)~W^)W>^ {p4r) ~ W^) W, 



(2.35) 



(2.36) 



where we have used the fact that p*(r) is an increasing function of r. Moreover at the boundary of the 
interval W{rQ — 5) = W {r^ + 5) = y^p^(l) which is not smaller than g^, thanks to the upper bovmd ( 2.26 ) , 
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which reads g^ < p^(l). Therefore W{r) is a supersohition to (2.31 1 in the interval [vq — S,ro + S] and by 
the maximum principle 



9*{ro) < W{ro) - y^ p.{ro + S) coth 
where we have used the fact that coth(a;) is a non-increasing function 



— ^2p^{ro + S) 



By the explicit expression (2.32) and the inequality (2.34), one has 

\p,{ro + S)-p4ro)\<Ce^n^S, 



h{r,) > P^"(ro) (l - Oil logel-i/^) > Cs^n^ir^ - Rl), 



so that (2.37) becomes 



g*{rQ)< Vp*iro) 1 + 



CS 



Hi 



coth 



3e 



^2p^{n, + 6) 



When the argument of coth tends to oo, i.e., 



-\/p^(?-o + '5)>l, 



we can bound 



and obtain the inequality 



coth(a:) = 



1 + e- 



1-e- 



^<(1 + Ce— ), 



ff.(^o) < ^hiro) ( 1 + 



C5 



i? 



1 + C exp 



25_ 
3e 



^J2p^{rn + 5) 



By (2.39) the second error term on the r.h.s. of the above expression is bounded from above by 



exp -C<5f7Jr2-i?2 , 



so that by taking 



6^Ce^\\oge\''{rl-Rl 



32\-l/2 



(2.37) 



(2.38) 
(2.39) 

(2.40) 



(2.41) 



(2.42) 



such an error can be made smaller than any power of e, since one can choose the constant coefficient in 
5 arbitrarily large. With such a choice the other error term becomes 

Ce^llogep {li-Rl)-"' > Ce^/^l loge^/^ 



since by definition r^ — i?^ < 1 — i?^ < ©(ej loge|). Hence the second factor on the r.h.s. of (2.41 1 can 

be absorbed in the above remainder. 

Moreover for any r > R^, S < e^^^j logej^^^ <C e^^^ and one can extend the estimate to any r > 

i?h + 0{e^^^\ logep). For the same reason the estimate applies also to the region [1 — 2(5, 1]: There one 

can use ( |2.26 l and the fact that /?*(1) — Pi,{l — 26) < CQ^e^S < Ce^^^, which is much smaller than the 

error term above. 

The final estimate is then 



gAr)<^M^)\l + Ce^\\ogemr'~Rl) '^'' 



(2.43) 
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for any R^ + Ce^/^l logep < r < 1. 

The next step is the replacement of p^, with p^^ by means of (2.33), which yields an additional remainder 
given by 



TF 



p- (r)-^ Mr) P'' {r)\ < C{en)-' (r^ - K) 



-1 II _ 



„TF| 



lL°=(^TF) 



< 



C e' \\og eY ir' ~ Ri 



(2.44) 



which can however be absorbed in the error term in (2.40), since r^ — i?^ < e| logej. 

In order to prove a corresponding lower bound, we fix some tq and 5 as before, i.e., such that 
i?h + Ce^'^^l logep +5<ro<l — 5 and (5 <C 1. Since p^ is an increasing function of r and g* is 
positive 

- Ag, > I [p,(ro - (5) - gl] g,. (2.45) 

Then we denote by h{r) the function solving for r £ B 

-Ah = r^ (1 - /i^) /i, 
with Dirichlet boundary condition /i(l) = and e — > 0. In [Se] it was proven that h satisfies the bound 

1 - cexp \ ^^ \ < h{r) < 1. 



2e 



If we now set 



h{r):=Vp4ro-5)h(^-^-^^ 



e := 



then h solves in [rp — (5, rg + S] the equation 



-Ah 



S^2Mro-Sy 



P*{ro -5) -h"^ 



(2.46) 

(2.47) 



with Dirichlet conditions at the boundary r — r()±5. Thanks to (2.45), 5* is a supersolution for the same 



problem, so that by the maximum principle gi,{r) > h{r) inside the interval and in particular 



9*{rQ) > h{ro) > \/p*(ro-(5) 



1 — c exp < 

^ 2e 



(2.48) 



for any i?h + Ce^/'^\ log ep + S < rp < 1-6. 



Note that by choosing 6 as in (2.42) the remainder in the above expression can be made smaller than 



any power of e. However the estimate of p*(r — S) in terms of p^(r) provides the same remainder as in 
the upper bound proof, i.e., 



g 



*{r) > \/h(r) 



l-Ce^\\ose?{r^-R 



ir'^ 



The extension of the estimate to the whole region [i?h + 0{e'^/'^\ logep), 1] as well as the replacement of 
Pi, with p^^ can be done exactly as in the upper bound and the remainders included in the above error 
term. D 



We conclude the section with an useful estimate of the gradient of the densities 5*: 
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Proposition 2.7 (Gradient estimate for 5*). 

^s e — > and for any i? > i?h + 0{e^/'^\ logep), one has 

l|V5.ILo=(6\8,) < ^£^' {R' - Riy'^' ll5*llL^(e) 

which inside A = {f : R^ + e| loge|^^ <''<!} becomes 

\\yg.h^^j)<Ce-y^\\oge\-'/\ 



Proof. We first exploit the fact that any g^, is radial to rewrite (2.22) as 

- 9': - r-'g: + m - u:fr-^g. + 2e-^ gl = (/i^P + 2U[n - u:]) g. 



(2.49) 



(2.50) 



(2.51) 



and take the L°° norm inside B \ Br, for some i? > i?h + 0{e'^/'^\ logep), of both sides to obtain (see 

<C\\r-'g:\ 



(2.32)) 



\9*\\L°°{B\Bn) 



\L°°{B\Br) 



+ Ce~^ Up,~gt)gJ 



< 



C 



L°-(B\Br 

2 (u-2 



2N-1/2 



9'.\\L-iB\BR)+^'^R'-Rl) \\9A\L-iB) ' (2-52) 



by (2.33) and the pointwise estimate (2.30): 

32\-l/2 



C (i?2 _ Rl)-'^' + 0(£-l/2| log£|-l/2) < C (i?2 _ Rl 



2\-l/2 



since i?2_ij2 < c'ellogel. 



On the other hand by the Gagliardo-Nircnberg inequality (see, e.g., [HI Theorem at p. 125]) 



1/2 



1/2 



II ' II <r' r"" II "II ' II II ' 

\\9*\\l'=°(b\Br) - ^ \\9*\\l'=°{b\Br) \\9*\\l'=°(b) 



< 



C 



32N-1/4 



1/2 



s~'{R'-Riy"^\\g4L^iB) + \\9'JL^iB\BR)\\9*\\L^iB) ' (2-53) 



which implies the result. 



D 



3 Reduction to an Auxiliary Problem on an Annulus 

The first step towards the proof of the main results is the reduction of the original GP energy functional 
to an analogous functional on a suitable annulus. The main ingredients of such a reduction are, on the one 



hand, the exponential smallness of the GP minimizers proven in the last section (see, e.g., (2.10), (2.11), 
(2.28), etc.), which intuitively implies that all the L^ mass and therefore the energy are concentrated in a 



annulus close to the boundary, and on the other a decoupling of the GP energy functional, which allows 



the extraction of the giant vortex energies introduced in (1.16). 



The second part of the section is devoted to the discussion of the optimal giant vortex phase: Whereas 
both the energy splitting and upper bound hold true for any reasonable phase w, they are useful only 
for specific choices of the phase. As we are going to see, the optimal phase luq can be defined as the 
minimizer of a suitable coupled problem in the annulus. 
We also discuss the existence of the analogous minimizer associated with the original problem in the ball 



B, i.e., the phase Wopt occurring in Theorems 1.2 and 1.3 as well as some relevant properties of it. 
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3.1 Energy Decoupling 

Before stating the main result of this section, we first recaU and introduce some notation. 
The main object of the reduction is an annulus 

A:={f : R<<r <1}, 



(3.1) 



where the inner radius i?< < i?h has to be chosen in a proper way so that two conditions are simultaneously 
fulfilled: The radius _R< should be sufficiently far from i?h in such a way that the exponential smallness 
proven in Propositions 2.2 applies inside the complement of A. At the same time, _R< must not be too 
far from i?h; in fact we shall in SectionHneed that |i?< — Rh\ <^ e|loge|~^ (see (4.89) and the use of 



(4.22) in the proof of Proposition 4.1 ) 



All these conditions are satisfied if one chooses 



R< :— Rh 



-e^l\ 



(3.2) 



The apparently strange power of e occurring in the above expression is essentially motivated by the 
pointwise estimates (2.111 and (2.291, which hold true for any r < i?h — ©(e^^^): In the definition of 



the inner radius of the annulus we could have chosen any i?< satisfying such a condition in addition to 
|J?< — R\\ <^ e\ loge|~^. In particular any remainder given by a power of e smaller than 7/6 but larger 
than 1 would have been all right. However for the sake of simplicity we make an explicit choice among 
the allowed remainders and pick £*/^. 

The auxiliary problem on the annulus is associated with the estimate of the energy functional 



£a,u,M •= 



A 



<irgXu;[\ 



Vvr - 2B^ 



(w,Vw)+e ^gX^ 



1- 



i-n'}, 



(3.3) 



where B^^ is defined in ( 1.32 1. According to the convention used in the rest of the paper, f^^ without the 



label A stands for the same energy as above but with the integral extended to the whole of B. 



A key ingredient in the proof of the GP energy asymptotics is a lower estimate of (3.3). In fact in the 



next Proposition, which is the main result proven in this section, we show that the energy (3.3) provides 
a lower bound to E'^^ — -E^^ when v is suitably linked to $^^: 



Proposition 3.1 (Reduction to an annulus). 

For any uj ^ Z such that \uj\ < 0{e^^) and for e sufficiently small 



^A,uj 



SaA^^] - 0{e°^) < E'^^ < E% + 0(e°°), 
where the function u^ £ H^{A) is given by the decomposition 



"^"-A^ -: g.A,c.(r•)M„(r)exp{^([^!] -c^)z?}. 



(3.4) 



(3.5) 



Proof. The proof of (3.4) is done by proving suitable upper and lower bounds to the GP ground state 
energy. 

The upper bound is obtained by testing S^^ on a trial function of the form g(r) exp{i([r2] — ijj)'d}, 



where g is an appropriate regularization of gA,ijj, and using the definition ( 2.19| ). Since gj,^^ does not 
vanish at i?<, it is not in the minimization domain of 5^^ and one has to regularize it at the boundary 
9i?<, e.g., taking 



gA^f) 

g(r):-c<( e-"g^,^(i?<)(i?< 




if r<£ A, 

if i?< -e" <r < i?<, 

if r<R<- £", 



(3.6) 
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where c is a normalization constant and n some arbitrary power greater than 0. Thanks to the exponential 
smallness ( 2.29[ ), which implies gA.uj{R<) = 0{e°°), the normalization constant is c = 1 — 0{e°°) and the 
energy of g in [i?< — e", i?<] is exponentially small as well. The upper bound trivially follows. 

The lower bound is mainly a consequence of a classical result of energy decoupling, which has been 
already used in different contexts, e.g., in [LM| . 
The starting point is however a reduction to the annulus A of the GP energy functional: Exploiting the 



exponential smallness (2.11) of ^^ , it is very easy to show that 



^GP ^ ^GP [^GP] > ^GP ^^GP] „ 0(goo)^ 



(3.7) 



where we have used the symbol £^^ to denote the restriction of the integration to A. Note that in the 
above inequality one can avoid the estimate of the gradient of \I>^^ by simply rewriting the functional as 



in (1.2) and neglecting all the positive terms given by the integration over B\A; the only negative term 



(centrifugal energy) can then be estimated by means of (2.11). 

Thanks to the exponential smallness, one also has that the L"^ mass of vJ/GP outside A is very small, i.e. 



* 



GPI 



\L^A) 



> 1-0(2°°). 



(3.8) 



Now inside A the decomposition (3.51 is well defined and therefore one can calculate 



2nm-uj)gXJuJ^ + e-'gX^{2Kf-l)} 
A simple integration by parts then yields 



since the boundary terms vanish because of the Neumann conditions satisfied by gA,uj at the boundaries 
OBr^ and dB. Then one can replace in the above expression AgA,uj by means of the variational equation 



(2.22) and the result is 



£T[^'''']=£AAu.]+f^T:J drg\ 

J A 



|2 -2 



I '^rgX^ 
J A 



-2 II ii4 



ii4(^) > SaA^u] + EX. - 0(£°°)> (3.9) 
by (3.8) and the definition of the chemical potential AS^ i^^^ Proposition 2.3 1. D 



3.2 Optimal Phases and Densities 

The idea behind the decomposition ( |3.5| is that, after the extraction from ^1^*^^ of a density g^^ and a giant 
vortex phase, i.e., the phase factor exp{i([r2] — uj)'d}, what is left is a function u^, which contains all the 
remaining vorticity in '^ . Therefore in the giant vortex regime, one would like to prove that |U(^| ~ 1 
inside a suitable annulus at the boundary of the trap, and the key tool to proving such a behavior is a 
detailed analysis of the reduced energy £a,u- 

However, in order to prove such a result, both the phase and the associated density gA,u have to be 
chosen in an appropriate way: The result stated in Proposition |3.1| is basically independent of w, i.e., it 
applies to any reasonable giant vortex phase w. Nevertheless the leading order term in the GP energy 



asymptotics is given in (3.4) by E^u which depends in crucial way on lo and it is clear that in order 
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to extract some delicate information about Ui^ like the absence of zeros, the estimate (upper bound) of 
the reduced energy £a.uj[uuj] through (3.4) has to be very precise. This leads to the definition of a giant 
vortex optimal phase (associated with the annulus A), which is nothing but the minimizcr of the energy 
-^A^uj with respect to ui. 

In the next Proposition we show that there exists at least one ujq minimizing E^^^ as well as some 
properties which are going to be crucial in the rest of the proof. 

Proposition 3.2 (Properties of the optimal phase uiq and density gA.i^o)- 
For every e > there exists an ujq ^ Z minimizing E^^^ ■ ^^ satisfies 

2 



Wo = 



S^/tte 



[l + 0{\loge\-'/')) 



(3.10) 



Moreover the minimizer gA,Luo '^f ^Aui satisfies the hound 

I <^rg\ 
J A 



f]-ti^^)^0(l). 



(3.11) 



Proof. The existence of a minimizing wo G Z can be easily proven by noticing that \mi^^±oo E'?^ = 



since by (2.19) 



^GP 



which implies that, for given e > 0, only finitely many oj G Z can minimize the energy. 



The main ingredient for the proof of (3.10) is the energy bound (2.21), which implies, for any a; G Z 
such that |a;| < Ce~^, 

E^: < eZ„ < eZ < eY + 0(e-^| logel-i), (3.12) 



by definition of ujq. Choosing now cj = [2/(3v^e)] in the r.h.s. of the above expression and using (A. 10 1 
for E'i^^ , we get the inequality 



E 



iTF 



LOo 



2 

3-v/7r£ 



^-2fi([r!] -O) < ^TF < E^"" + ^^2nm ~ n) + Ce-^\\oge\-' (3.13) 



which yields the result. 



We now prove (3.11). A simple estimate yields 



9Aa 



£aIo±i [9a.uo] = ^aIo [9a,uo] ±2 f df[n- m - coo) r-'] g%^„ + j dr . 

J A J A 

EX.o ±2 f dr[n- m - wo) r-'] gX^o + ^<'- (3-14) 
J A 

Now suppose that 



A 



> 



3' 



df[n-i[Q]-uo)r-^]gX^, 
then, since i?< = 1 — o(l). ( |3.14 | would imply that there exists some a; G Z equal to wo ± 1 such that 

pGP ^ cGP r 1 ^ pGP 

^A^oj ^ '-'Am [9A,ujo\ ^ ^A,ujo' 

which contradicts the fact that ujq minimizes E^^^- Note that the above argument implies that the 



constant on the r.h.s. of (|3.1l|) is actually smaller than 1, i.e., 

dr[n-i[n]-ujo)r-^]gX. 



A 



< 1. 



(3.15) 

D 
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As anticipated above and in the Introduction, there is another giant vortex phase that naturally 
emerges in the study of the above functionals, i.e., the one associated with the energy in the whole ball 
B: More precisely it is defined as the minimizer of 

i^r:= inf e^i/] 

1/ II2— 1 

with respect to w € Z and denoted by LJopt'- 

Proposition 3.3 (Optimal phase ujopt)- 

For every e > there exists an Wopt G ^ minimizing E^^ . It satisfies 



t^opt 



Sv^Tre 



(l + 0{\\oge\-'/')). (3.16) 



Proof. The proof is basically identical to the one of (3.10): For instance it is sufficient to replace (3.121 



with the corresponding version in B (see (2.21 1). D 



4 Estimates of the Reduced Energy 

In this section we study the auxiliary problem introduced in Section l3] From now on we shall simplify 
notation by dropping some subscripts: 

QA.uo ■= 9, 
also 

B(r) := B^,ir) ^ ff^r - It^) e",, (4.1) 



and 

r,4 



£[v] := / dfL^\Vv\^ ~2g^B-{iv,Vv) + ^{l-\v\^] 
The following energy is crucial in our analysis: 

^H:=|^dr|5^|V.r + fJ(l-|«n^}. 

We recall that u := u^^g is defined by 

*GP(-) ^. g(r)u{r)exp{i{[n]~ujQ)'d} (4.3) 

and that from Proposition |3.1| we have 



(4.2) 



£[u] < 0{e°°). (4.4) 

We also need to define a reduced annulus 

A:^{f : i?> < r < 1} (4.5) 

with 

R>:=Ri, + e\logs\-\ (4.6) 



An important point is that from (|2.30|) we have the lower bound 

C 

ellogeP 

The main result of this section is 



g\r)>— ^onA (4.7) 
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Proposition 4.1 ( Bou nds for the reduced energies). 

Let u be defined by {4-3). If flo > 2(37r)^^ we have for e small enough 



T\u] < C 



|loge|5/2 



£[u] > -C — 



ei/2(log|loge|)2 
loge|3/2 



£1/2 log I logej 



(4.8) 
(4.9) 



Since the proof of Proposition |4.1| is rather involved we sketch the main ideas before going into the 
details. 
It is useful to introduce a potential function F defined as follows: 



F{r):=2 ds g^ {s) ( ^s - {[n] - ujq) 



l)-< 



ds g^{s)B{s) ■ e^. 



We have 



V-^F =: 2g'^B, 



F{R<) = 0, 



(4.10) 



(4.11) 



i.e., F is the "primitive" of 2g^B vanishing at i?<. We refer to Subsection 4.1 for further properties of 
F. Integrating by parts we have 

/ df\g^\Wu\^-2g'^B-{iu,Wu)]= / df\g^\Wu\^+Fcm\{iu,Wu)]- / da F{l){iu,dru) (4.12) 
Ja ^ ^ Ja ^ ' JdB 

and thus the energy £[u\ can be rewritten as follows 

£[u]= I dfg2|Vu|V / df Fcurl(TO, Vw) - / d<7 F{X){iu,dru) + [ df^{l-\u\^f 

J A J A JdB J A £ 



(4.13) 



It is in this form that the energy is best bounded from below. 



The boundary term (third term in (4.13)) is estimated in the following way: The property of Wq given 
in (3.11 ) implies that F{\) = 0(1). We combine this fact with an estimate of the circulation of u on the 



boundary of the unit ball that we provide in Subsection 4.2 Proving this estimate requires to derive a 
PDE satisfied by u and use it in much the same way as in the proof of the Pohozaev identity [P] . 
The first two terms can be estimated in terms of the vorticity of u. Indeed, suppose that |u| ~ 1 except 
in some balls (that we identify with vortices) whose radii are much smaller than the width of A. Let us 
denote these balls by {B{aj,t)} . j with J C N and t <C e| logej, i.e., much smaller than the width of the 
annulus. Then by Stokes theorem, if the degree of u around dj is dj (we systematically neglect remainder 
terms in this sketch) 

/ df Fcurl(iM, Vu) ~ V 2TrF{aj)dj. (4.14) 

Minimizing the sum of the first and the last term with respect to t yields t ex e'^/2| logej 1/2 (see also |CY[ 
p. 6] for heuristics about the optimal size of the vortex core), which implies an estimate of the form 



/ drg'\Vuf>Y,2ng^ia,)\d,\\og('^^)>J2^9'ia,)\d,\\log 



(4.15) 



In Subsections |4.3| and |4.4| we give a rigorous version of this heuristic analysis. The main tools were orig- 
inally introduced in the context of Ginzburg-Landau (GL) theory (we refer to |BBH21 [S5] and references 
therein). The method of growth and merging of vortex balls introduced independently by Sandier |Sa| 
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and Jerrard [JT] provides a lower bound of the form (4.151 (see Subsection 4.3 1. On the other hand, the 
jacobian estimate of Jerrard and Soner |JS| allow to deduce in Subsection 4.4 that 



curl(w, Vu) ~ 2. 2ndjS{'r~ Qj), 



(4.16) 
i.e., the vorticity measure of u is close to a sum of Dirac masses. This implies an estimate of the form 



(4.14). 



Having performed this analysis, the role of the critical velocity becomes clear: We have essentially (note 
that the boundary term is negligible in a first approach) 



£[u]>J227r\d, 



-y{a,)\\oge\+Fiaj) 



(4.17) 



If flo > 2(37r)~"'^ the sum in the parenthesis is positive for any Sj in the bulk (see the Appendix), which 
means that vortices become energetically unfavorable. 

So far we have assumed that the zeros of u were isolated in small 'vortex' balls. One can show that this 
is indeed the case by exploiting upper bounds to J^[u] (this energy controls, via the co-area formula the 
size of the set where \u\ is no close to 1). We derive a first bound from our a priori estimate on £[u] 
(4.4) in Subsection 4.1 We emphasize however that this first bound is not strong enough to construct 



vortex balls in the whole annulus A. To get around this point we split in Subsection |4.3| the annulus into 
cells and distinguish between two type of cells. In 'good' cells we have the proper control and perform 
locally the vortex balls construction. On the other hand we are able to show that there are relatively 
few 'bad' cells where the construction is not possible. Thus the analysis in the good cells allows to get 
lower bounds for £[u]. These bounds can in turn be used to improve our control on J"[u] and reduce the 
number of bad cells. The analysis can then be repeated, but now on a larger set. Finally we see that 
there are no bad cells. Such an induction process is 'hidden' at the end of the proof of Proposition 4.1 in 



Subsection 4.5 (see in particular the discussion after (4.101)) 



4.1 Preliminaries 

We first give some elementary estimates on B and F that we need in our analysis: 
Lemma 4.1 (Useful properties of B and F). 



Let B and F be defined in (4.1) and (4.10) respectively. We have 



\B\ 



L-^{A) 



-^IIl~(^) - ^^ ' 



IVFI 



L^[A) 



\F{1)\^2 



ds g^{s)B{s) ■ e^ 



< Ce-^\loger\ 



< C. 



Moreover there is a constant C such that 



|F(r)|<Cmin' l*^ ■^<'"2 



(r), 1 + 



C 



e2|loge| 



r-l 



(4.18) 
(4.19) 
(4.20) 

(4.21) 



(4.22) 



for any f d A. 
Proof. We have 



B{r) = {{n - [n])r + [n]{r - r-i) + t^or"^) e^ 
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so (4.18) is a consequence of I Wo I < Ce ^ (see (3.10)) and the fact that |^| (xe|loge| (see (3.1) and (3.2)) 
and i} oc e~^\ log£|~^ 



We obtain (4.20) from (2.25), (4.11) and (4.18). Onthe other hand F(i?<) ^ 0, which imphes that (4.19) 



follows from (4.20) and |.4| oc e| logej. Moreover (4.21) is exactly the same as (3.11) in Proposition 3.2 



To prove (4.22) we prove that |-F'(7')| is smaher than both terms on the right-hand side. From equation 



(4.10) we have 



|F(r)| < 2||i?||^^(_4) r dsg'is) < Ce-' f ds g' 



is) 



but g^ is increasing so 



As g\s) < g\T)\r - R<\. 



We deduce that |-F'(r)| is smaller than the first term on the right-hand side of (4.22). On the other hand, 
using (4.20) and (4.21) we have immediately 



\F{r)\<C{l + e-^\\oge\-\-l\) 



for some finite constant C. Thus (4.22) is proven 



D 



We now prove the energy bounds that are the starting point for the vortex balls construction in 
Subsection W. 



Lemma 4.2 (Preliminary energy bounds). 

We have, for e > small enough, 



F[u] < 



C 



£\u] > --S. 



Proof. We have, using (4.18) together with the normalization of ^'^^, 
2 



dfg B ■ {iu,Vu) 



1 



< - / drg^ |Vur + 2 / df 5^B^|u|^ < - / df g-" iVuf + 



C 



(4.23) 

(4.24) 

(4.25) 



Equation (4.24) immediately follows. We obtain (4.23) from (4.4): 



T[u] = £[u] + 2 / dfg-'B ■ {iu, Vu) < 0{e°°) + 2 
Ja 



dr g B ■ {iu, Vu) 



< 



i^drV|Vu|V^<^-FM + ^. (4.26) 



Subtracting ^^^[m] from both sides of the last inequality we get (4.23) 



n 



4.2 Equation for u and Boundary Estimate 

We first derive from the equations for ^'^^ and g an equation satisfied by u: 

Lemma 4.3 (Equation for u). 

Let u he defined by (4.3). We have on A 



V{g^Wu) - 2ig^B ■Wu + 2^ {\u\^ - l) u = Xg'^u, 



(4.27) 
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where A G M satisfies the estimate 

Moreover u satisfies the boundary condition 

drU — on dB. 
Proof. In this proof we use the short hand notation 



The equation for u is a consequence of the variational equation (1.5 1 satisfied by '^'^^ 



A*GP - 2n ■ L*GP + 2£-2 I^GPI- ^GP ^ ^GP^ 



and that solved by g (see (2.22)), 



VL 



I\g+\nr ) g- Q^r^ + 2e"^g-' = /i 5- 



Direct computations starting from (4.3) show that 

,GP 



-A*^ 



n^ 



and 



uAg — gAu — 2Vm • Vg H — :vgu -\ — ^g^ ■ Lu I e 



2n ■ L-^^^ =i- 2gn ■ Lu - 2nClgu\ e'^^ . 



iQ.-d 



Plugging the equation for g in (4.331 we obtain 

- A*^P = f-gAu - 2Vu • Vg + jj!^^ gu + 2Q.Q.gu - 2e^^g^u + 2gCl ■ Lu\ e' 



n§ 



Next, combining (4.311, (4.34) and (4.35), we have 



where 



2 2 ^ ^ ^ ^ 

g/S.u - 2Vu • V5 + —g^ (luP - l) u + — yfi • iu - 2,gf7 • Lu = Agu 



GP r.GP 



A = /i^^ _ /t' 



There only remains to multiply (4.36) by g and reorganize the terms to obtain (4.27) 



(4.28) 
(4.29) 
(4.30) 
(4.31) 

(4.32) 

(4.33) 
(4.34) 

(4.35) 

(4.36) 

(4.37) 



The Neumann condition ( 4.29 ) is a straightforward consequence of the corresponding boundary conditions 
for g and ^^^ . 



(3.8), 



We now turn to the proof of (|4.28|. Multiplying (|4.27|) by u* and integrating over A, we obtain, recalling 

(4.38) 



|A|(l-0(e-))<|A| / df|v|/GP|2^|A| I drg^\u\'^£[u]+J df'^{\u\*-l) 



A 



.9 



Using Cauchy-Schwarz inequality 



A <^ 



< 



(/.?(I"I^-'>1"'(/.?<'''I"I^ + '''' 



1/2 



so (4.28) follows using the upper bounds (2.2) and (2.25) on 5 |u| — \^!^^Y ^^d g respectively together 
with \A\ cxe|loge|. D 
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We remark that the equation satisfied by u is exactly of the form that we would have obtained if u 
had been a minimizer of £ under a mass constraint for gu. 

We are now able to estimate the circulation of u at the boundary of the unit ball: 

Lemma 4.4 (Boundary estimate). 

We have, for £ > small enough, 



/^^dag^|aHV/^^d.g;(l-H^)^< ,3/.|w.|i/. -^M + g|gM 



£3/2 1 log £ 1 1/2- ^^J ' e'^^-JI- (^-^^^ 

Proof. The main idea is a Pohozaev-like trick, like in |BBH2[ Theorem 3.2]. Wc multiply both sides of 
(4.27) by (r — _R>) eV • Vu*, then integrate over A C A (see (4.5) for the definition of A) and take the 
real part. Recall that 1 — i?> oc e\ log£|. 
The first term is given by 

-- / df{{r~Ry)er-Vu*V{g^Vu) + ir-Ry)er-VuV{g^Vu*)} = 

drh (l - ^) \Vuf +g'^ |e, • Vuf + \g^ {r - i?>) e, • V|Vw|2 



A 

If dag^l- i?>) \dru\' - f dfig [r - i?>) e. • Vg |Vwp ' ^> -^ 

^ JdB J A L 



2r 



5 



2|e;- VmP-IVwP 



(4.40) 



where we have integrated by parts twice, remembering the Neumann boundary conditions for g and u on 

dB. 

The third terms yields 

\ I df{(r-i?>)e;- Vu*/(|up-l)u+(r-i?>)e;- Vw/(|up-l)w*} = 

-i^/drV(r-i?>)e.-V(l-H2)' = -l^ /■ da (1 - i?>) .g^l - l«P)' " 

^^df|(^l-|^)/(l-H2)V25='(r-i?>)e.-Vg(l-H2)'|. (4.41) 
Altogether we thus obtain 

\{l-R>) f da(5^|9H' + 4(l-|"l')' 



df<^5 (r-ii'>)e:, -V^IVmI- 



i? 



2r 



> „2 



2|e;. -Vul^- IVmI^ 



2i.g^ (r - i?>) eV ■ Vu* B ■ Vu 



^J^drl^(^l-^^g^l-\u\^y + 2g'{r-R^)e,-Wg{l-\u\^fy 

^jdrg\r-Ry)er-V\u\\ (4.42) 



We then estimate the moduli of the terms of the r.h.s.. Obviously 

^4 / R^ 



1^] ^> 2 



2|e;- Vup-lVwP 



£^ 



1- 



2r 



(i-H^)' <-FM- 



Using (4.6) and (4.18) we obtain 



drig"^ (r - i?>) eV • Vm* B • Vu 



<C|log£|/ dfg-'lVur <C|l0g£|J'[M] 



(4.43) 



(4.44) 
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whereas, using the normaUzation of ^"^^j (4.23) and (4.28), 

/ df.g^(r-i?>)e; • V|u 
Ja 



A / dfg u[r ~ i?>) eV • Vu* 



^^Uog<I^MI + ^37^^^M^/^)(//' 



r.g2|„|2 



< 

1/2 



1/2 



C|log£||£:[u 



C|l0g£|l/2 



rl/2 



J"[u]. (4.45) 



Combimng (2.50) with (4.7), we obtain 



^ j drc/{r~Ry)er-Vg (l 
J A 



< 

C|l0g£|'^/4 

^174 



^ 



dr^9\l^\u\r< 



C|l0g£|9/4 



^1/4 



By similar arguments 



/ dfg (r-i?>)e^-V5|Vu|- 



< 3^^H^ 



-F^. (4.46) 



(4.47) 



1 - \u\^f [ < C\ logel \£[u]\ + ^'^°f^f V m 



Gathering equations (4.42) to (4.47), we have 

hl-R>) f daig'\dru\' + ^4^ 
2 JdB I £ 

and there only remains to divide by 1 — Ry ex e\ loge| to get the result. D 

4.3 Cell Decomposition and Vortex Ball Construction 

In this subsection we aim at constructing vortex balls for u. Namely, we want to construct a collection of 
balls whose radii are much smaller than the width of A and that cover the set where ImI is not close to 1. 



The usual of way of performing this task is to exploit bounds on T[u]. Unfortunately, the bound (4.23) 
is not sufficient for our purpose. It only implies that the area of the set where u can possibly vanish is of 
order e^| logep, whereas the vortex balls method requires to cover it by balls of radii much smaller than 



the width of A, which is 0(e|loge|). However, the estimates of Lemma 4.2 are definitely not optimal 



and can be improved by using a procedure of local vortex balls construction. 



The idea is the following: If the bound (4.23) could be localized (in a sense made clear below), we could 
construct vortex balls. As this is not the case we split the annulus A into regions where the bound can be 
localized and therefore vortex balls can be constructed as usual, and regions where this is not the case. 

Definition 4.1 (Good and bad cells). 

We cover A with (almost rectangular) cells of side length Ce|loge|, using a corresponding division of 
the angular variable. We denote by N (x e^^|loge|^^ the total number of cells and label the cells as 
An,n G {1, ..., A''}. Let < a < ^ be a parameter to be fixed later on. 



We say that An is an a- good cell if 



J^ df|5^|Vz.r + fJ(l-H2)^} 



< lioggL - 



(4.48) 



We denote by N^ the number of a-good cells and GSa the (good) set they cover. 



Giant Vortex - CRY- December 6th, 2010 



32 



We say that An is an a-bad cell if 






,2 1 ^ llogeL - 



(4.49) 



We denote by N^ the number of a-bad cells and BSa the (bad) set they cover. 



Note that the annulus A has a width ^ ex e| loge| (which imphes that N ex e^^\ loge|~^) so that we 
are actually dividing it into bad cells where there is much more energy that what would be expected from 
the localization of the bound (4.23) (namely Cie~^ ex e~^\ logej) and regions (good cells) of reasonably 
small energy. A first consequence of this is, neglecting the good cells and using ( 4.23[ ), 



N^< 



loge 



■e"^M < 



C 



£|log£ 



■£" < N, 



(4.50) 



i.e., there are very few a-bad cells. A consequence of the refined bound (4.8) that we are aiming at is 
that there are actually no a-bad cells at all. 

We now construct the vortex balls in the good set. The proof is merely sketched because it is an 
adaptation of well-established methods (see [SS] and references therein). Note that the construction is 
possible only in the subdomain A where the density is large enough. 

Proposition 4.2 (Vortex ball construction in the good set). 

Let < a < 2- There is a certain Eq so that, for e < Eq there exists a finite collection {Biji^j := 
{B{ai, gi)}j^^j of disjoint balls with centers di and radii Qi such that 

1. {reG5„ni : \\u\-l\>\\oge\-^](Z{j^^,B,, 

2. for any a-good cell An, Z]j,e,n^„#0 ^ ^ '^l logel"^- 

Setting di :— deg{u^dBi}, if Bi <Z AC] GSa, and di — otherwise, we have the lower bounds 



I drg' \Vu\' > 27r Q - a\ \d,\g'{a,) |log£| fl 



-C 



log I log £ I 
|loge| 



(4.51) 



Proof. We begin by covering the subset of GSa H A where u can possibly vanish by balls of much smaller 



radii than those announced in the Proposition. Let An be an a-good cell. We have, using (4.7) and 



(4.48), 



A„nA 



dfi iVul 



£^1 log£P 



< 



C£|log£|3 f dfL^\Vu\ 



"^^^ <C|l0g£|4£-". 



Then, using the Cauchy-Schwarz inequality, 

^jvkiili 

dr 



lA„nA £3/2|l0g£|3/2 

and the coarea formula implies 

|l-f2' 



<C|l0g£|'*£-" 



tei 



+ '^^£3/2|l0g£|3/2 



n^(^{reAnnA : |M|(f)=t}) <C|log£|V", 
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Figure 1 : Initial collection of vortex balls 



where Ti^ stands for the one-dimensional HausdorfF measure. We argue as in |SS| Propositions 4.4 and 
4.8] to deduce from this that the set {r G An^iA : ||u| — 1| > | loge|~-^} can be covered by a finite number 
of disjoint balls {B{aj, Qj)}j^j with X],6j 9j — Ce^^^~"\ loge]^'^/^. Doing hkewise on each a-good cell, 
we get a collection {S(aj, Qi)}iei, covering {re GSa ri A : ||w| — 1| > | loge|^^} and satisfying 



E 



Qt 



< C£3/2-"| log£|l^/^ for any a-good cell A„ 



(4.52) 



i<£l,B{ai,gi)nA„^0 



The balls in this collection may overlap because we have constructed them locally in the cells. However, 
by merging the balls that intersect as in [SS) Lemma 4.1], we can construct a finite collection of disjoint 
balls (still denoted by B{ai, Qi)) satisfying the same bounds on their radii and still covering {r € GSaHA : 

||w|-l|>|l0g£|-l}. 

The collection B{di, Qi) is represented in Figure 1. The gray regions are those were u could vanish, i.e., 
bad cells and vortex balls. The dashed circle is the inner boundary of A. 

We now let the balls in our initial collection grow and merge using the method described in [SS| Section 
4.2], adding lower bounds on the conformal annuli constructed. In brief, the idea is to introduce a process 
parametrized by a variable t (interpreted as 'time') so that, as t increases, the balls will grow all with the 
same dilation factor. There are two phases in the process: When balls grow independently their radii get 
multiplied by some factor, say m{t), and we add lower bounds to the kinetic energy of u on the annuli 
between the initial and grown balls. The first time two (or more) balls touch we merge them into larger 
balls (see ^SS. Lemma 4.1]), and continue with the growing phase. Figure 2 shows this process. 
The lower bounds on the annuli are obtained by integrating the kinetic energy over circles centered at 
the balls centers Si during the growth phases: 



/ 



BiSi,m{t)e^)\B(ai,e^) 



df |Vup > 271^2 log (m(i)) (l - | loger^) 



(4.53) 



where di := deg{u, dB{ai, Qi)}. 

The main point in this computation is that ||m| — 1| < | loge|^^ (in particular u cannot vanish) between 
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Figure 2: Growth and merging process 



the circles dB{di, Qi) and dB{ai, m{t)gi). Thus the degree of u is di on any circle dB{di, g) with gi < g < 

m{t)g^. 

To add the lower bounds obtained during the different growth phases we use that if two balls, say B{dp, gp) 

and B{dq, gq) merge into a new ball S(as, gi^) one always has dp + d^ > |dp| + |(iq| > |(ip + (ig| = |ds|- This 

accounts for the fact that we get a factor \di\ in (4.51 ), whereas the factor in ( 4.53| ) is df . 

We stop the process when we have a collection of disjoint balls {Bi}if=i := {6(0^,^)},-^^ satisfying 
condition 2 of Proposition |4.2| (property 1 is satisfied by construction, since it holds true for the initial 
family of balls), i.e., the balls are large enough, so that we have 



BiSi.Qi 



dr\Vuf>2n{l^a)\d.\\lo,e\(l-c'^^ 



(4.54) 



where di := deg{u,dBi}, if Bi C GSa H A^ and di = otherwise. The final configuration is drawn in 
Figure 3. Note that the radii of the final vortex balls are still much smaller than the diameter of the cells. 
This is important for ( |4.55[ ) below and the jacobian estimate in the next section. 

The logarithmic factor (5 — a) |loge| ( 1 — C "mq ^| ) comes from the logarithm of the dilation factor of 

the collections of balls, i.e., 






> C£"-1/2M |-23/2^ 



On the other hand, using (2.301 and the fact that |Vp \ <Ce \ loge| ^ (see (1.11)), we have 



m\ng^{r)-g^{ai) 

rGBi 



<C £i/2|loge|3/2||pTF| 



L=°iA) 



minp (r) - p (a^) 
f&Bi 



< 



C (e-i/2| iogg|i/2 ^ e-2| loge|-2g,) < Ce-'\loge\-' < C\ logerV( 



(4.55) 



because di & A and Qi < Cejlogej ^. We conclude from (4.541 and (4.55) that the lower bound ( |4.51 1 
holds on each ball we have constructed by bounding below g'' with its minimum on the ball Bi. The error 



minpgBi 5 (r) ~ 9 V^i) can then be absorbed into the g (ui) log | loge| term. 



D 
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Figure 3: Final configuration of cells and vortex balls. 



4.4 Jacobian Estimate 

We now turn to the jacobian estimate. With the vortex balls that we have constructed, it is to be 
expected that in the a-good set the vorticity measure of u will be close to a sum of Dirac masses 

curl(iu, Vu) ~ 2. '^'n'di5{r — Ui), 

where S{f— Si) stands for the Dirac mass at Oj. Indeed, outside of the balls |u| ~ 1, so curl(m, Vm) ~ 0, 

and the balls have very small radii compared to the size of A. 

Proposition |4.3| gives a rigorous statement of this fact. Note that we can quantify the difference between 

curl(iu, Vu) and '^^ 2TTdiS(f— Si) in terms of the energy only in the domain A where the density is large 

enough. 

Proposition 4.3 (Jacobian estimate). 

Let < a < 2 o,nd cj) be any piecewise-C^ test function with compact support 

supp(0) C ^n GSa- 



Let {Sijjgj := {B{di, gi)}^^j be a disjoint collection of balls as in Proposition ^.2. Setting di 
deg{u,dBi}, if Bi C An GSa, and di = otherwise, one has 



> 2tt di(j){di) — I dr0curl(m, V 

i^i Jas^nA 



u) 



<C\m\\^^^^s. 



e'\\0Ke\-'T[ 



(4.56) 



Proof. We argue as in |SS[ Chapter 6]. We first introduce a function ^ : M+ -^ M+ as follows: £,{x) = 2x, 
if x e [0, 1/2], and £_{x) = 1, if a; € [1/2, +oo[. This function satisfies 

• ^{t) < 2t and ^'(i) < 2, 

• |^(i)-i|<|l-t|and|e(t)-l|<|l-t|. 



Giant Vortex - CRY- December 6th, 2010 



36 



• |eW"i'| <3t\l-t\, 

and we define w as a regularization of u (in tlie sense that |w| = 1 and therefore curl(iw, Vw) — when 
|u| is far enough from 0): 

w = u. (4-57j 

|u| 

Remark that we have (iw, Vw) — j^[iu^ Vw) and this has a meaning even if u vanishes. By integrating 
by parts on GSa H A and using the assumptions on (/>, one has 

/ df [curl(iw, Vu) — curl(ii(;, Vw)] = — / dr [(iu, Vu) — (ity, Vu;)] V^0. (4.58) 

JGS^nA JGS„r\A 



Now, 



/ dcr [(iu, Vu) — (iw, Vw)] V" 

JGS^nA 



< llW|lioe(Gs„) / dr|(zM,VM)-(zw;,Vw;)| < 

J A 



IV0I1 



L'^iGS^) 



|u|2-|u;|2 



|Vw| < \\y4>\\L^(Gs^) 111 - HWmA) W^Alha) < 



CIlVc 



i~(GS„) 



r5/2|loge|^/2^M (4.59) 



using the properties of ^, the Cauchy-Schwarz inequality, the definition ( 4.2 ) of T[u] , g^ > Ce ^ | log e| ^ 
on i and (1 - |w|)^ < {l - \u\'^f . 
We now evaluate 

/ dr curl(iu', Vt(;)(/) = > / dr curl(iw, Vw)0, (4.60) 

which follows from the fact that ||u| — 1| < | loge|~^ outside Ui^iBi, so that \w\ = 1 and curl(iw, Ww) = 
outside UifziBi. If Bt c An GSa, we have 



df |curl(iw,Vu;)||(/)(f)-0(a,)l <C||V(/.||i^(s.)ft / df|Vu|' 

B, JBi 



(4.61) 



and 



/ dfcurl(iw, Vw) = deg{u, ^B^} = 2ndi (4.62) 

by definition of w and d,;. On the other hand, if Bi ^ ^ n GS^ then Bid d (An GSa ) ¥" ^ ^^'^ thus 

/ dr curl(iw, Vw)0 < / df |(/)||Vu;|^ < 4 / df |0||Vm|^ < 

JB^nA JBinAnGS" JB.nAnGS" 

CIIV 



L^m^^l df|Vu|^ (4.63) 

'B.nAnGS^ 



(4.63), we obtain 



because |Vw| < 2|Vu| and 4> is supported in the interior of ^ n GSa- Gathering equations (4.58) to 

Ar (f) curl(m, Vw) 



y^(ij0(a,) 



iei 



GSa,r]A 



<C||V0|L^(c^^)£^/2|loge|9/2^ 

CY.\\^^\\L^i.B.)e^ 



Bi^AnGSc 



df|Vup (4.64) 
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and the result follows because, using g.i < e\ loge| ^ and (4.7) 



BinAnGS^ 



df |Vup < C||V0||^^(cs„)e'|log^r'E / drVlVup < 

2|l„„^|-2- 



C||V0||i^(GS^)e^|loger^-FM. (4.65) 

D 



Note that (4.561 is equivalent to saying that the normof curl(TO, Vu)— ^^ d,;(5(r— a^) in {Cl{Ar\GSa))* 



i.e., the dual space of C^{Ar] GSa), is controlled by the energy. 



4.5 Completion of the Proof of Proposition 4.1 



We now complete the proof of Proposition |4.1[ collecting the estimates of the preceding subsections. 
We want to avoid any unwanted boundary term when performing integrations by parts in the proof below. 
Indeed, our radial frontiers between the good set and the bad set are somewhat artificial and have no 
physical interpretation. Therefore it is difficult to estimate integrals on these boundaries. 
To get around this point we need to introduce an azimuthal partition of unity on the annulus in order to 
'smooth' the radial boundaries appearing in our construction. This requires new definitions: 

Definition 4.2 (Pleasant and unpleasant cells). 

Recall the covering of the annulus A by cells Am n £ {1, ..,-/V}. We say that An is 

• an a-pleasant cell if A„ and its two neighbors are good cells. We note PSa the union of all a- 
pleasant cells and N^ their number, 

• an a-unpleasant cell if either An is a bad cell, or An is a good cell but its two neighbors are bad 
cells. We note UPSa the union of all a-unpleasant cells and N^^ their number, 

• an a-average cell if An is a good cell but exactly one of its neighbors is not. We note ASa the union 
of all a-average cells and N^ their number. 



Remark that one obviously has, recalling (4.501 



" - 2 " 



<^N 



(4.66) 



and 



N^ < 2N^ < N. 



(4.67) 

The average cells will play the role of transition layers between the pleasant set, where we will use the 
tools of Subsections |4.3| and |4.4[ and the unpleasant set, where we have little information and therefore 



have to rely on more basic estimates (like those we used in the proof of Lemma 4.2 1. To make this precise 
we now introduce the azimuthal partition of unity we have announced. 

Let us label C/PS'„,/ e {1,...,L}, and PS^,m G {1,...,M}, the connected components of the a- 
unpleasant set and a-pleasant set respectively. We construct azimuthal positive functions, bounded 
independently of e, denoted by xf and Xm (the labels U and P stand for "pleasant set" and "unpleasant 
set") so that 

1 on UPSl, 

on PS*™, Vm e {1, 

1 on P5™, 



Xi 

xf 

p 



X'n 



. ,M}, andont/PS'i', W ^l, 



/ J Xm 



E 



Xi 



on UPSl, \fle{l,...,L}, and on PS*™', Vm' ^ m, 
= 1 on A. 



(4.68) 
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It is important to note that each function so defined varies from to 1 in an average cell. A crucial 
consequence of this is that we can take functions satisfying 



vxH < 



c 



e|loge| 



Vx„ 



< 



c 



el log el 



(4.69) 



because the side length of a cell is ex e| loge|. For example one can choose this partition of unity to be 
constituted of piecewise afRne functions of the angle. 
We will use the short-hand notation 



M 



Xin • 



Xout 



Am: 



E 

m— 1 



(4.70) 
(4.71) 



The subscripts 'in' and 'out' refer to 'in the pleasant set' and 'out of the pleasant set' respectively. 

We would like to use the jacobian estimate of Proposition |4.3| with cj) = XinF, whose support is not 
included in A but only in A (moreover it does not vanish on dB). We will thus need a radial partition of 
unity: We introduce two radii R^^t and R, 



Re 
r: 



as 

:= l-£|loge|-i, 
-cut - R>+e\\oge\~^ 
Let $in(?') and ^out(^) be two positive radial functions satisfying 

Cin(r) 

6„(0 

S,out{r) 

sin I Sout 



(4.72) 
(4.73) 



= 1 for i?-,t < r < i?+t, 

= for i?< <r< Ry and for ?' 

= 1 for i?< < r < i?>, 

= for i?,-^ < '^ < ^c^ut, 

= 1 on A. 



(4.74) 



For example ^in and ^out can be defined as piecewise afiine functions of the radius. Moreover, because of 

C|l0g£| 



(4.72) and (4.731, we can impose 

IVCin 



< 



|VCo 



< 



C|loge| 



The subscripts 'in' and 'out' refer to 'inside A' and 'outside of A' respectively. 

In the sequel {Bi}^^j :— {B{ai, Pi)}ig/ is a collection of disjoint balls as in Proposition 
sake of simplicity we label Bj, j ^ J C /, the balls such that Bj C .4 n GSa- 

Proof of Proposition \4-l\ 



4.2 



(4.75) 



For the 



Recall the properties of i^ (4.11 1. By integration by parts, we have 



I drL^\Wu\^ -2g'^B ■{iu,Wu)\ = f df^g^ \Wu\^ + Fcm\{iu,Wu)\ - f da F{l){iu,dru) (4.76) 
and we are going to evaluate the three terms using our previous results. We begin with a lower bound 



on the kinetic term in (4.76), using Proposition 4.2 We introduce a parameter 7 to be fixed later in the 
proof and estimate 

/ dfg2|Vu|2>(l-7)V / dfCi„g2|Vu|2 + (l-7) / drXutff'|Vu|2 + 7 / dfg^\Vu\\ (4.77) 

J A ,,- T JBi J A J A 
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Using the lower bound (4.511, we have 



JBj JBj V<^t}, ) Jg. 

'^°^'^°^^'^ ^^/d.V|V.p. (4.78) 



27r - - a |dj|6„(a,y (aj) |loge| 1 - C 



|loge| 



C 

|log£|- Je, 



The estimate of infg^ ^jn — Ciii(oj) is a consequence of (4.75) combined with Qj < Ce\ loge| ^. 
We now compute 

/ dfF curl(w, Vu) := / dr [^inXm^ curl(m, Vu) + CoutXin-F curl(TO, Vu) + Xout-F curl(iu, Vu)] . 

(4.79) 
We can use Proposition 4.3 to estimate the first term because ^inXin^ is a piecewise-C^ function with 



support included in ^ n GSa ■ We obtain 

/ df^inXinF cm\{iu,Wu) > 27r^ djF(aj)^in(aj)xi„(aj) - C ||V(^inXin-F)|li,oo(Gs^) e^| logep^J^M. 



j&J 



(4.80) 



Now, using (4.19), (4.20), (4.69) and (4.75), we have 



|V(6„XinF)|U.(.) < ^^. 



so that 



/ df^inXinF curl(m, Vw) > 27r V" djF(aj)^i„(aj)xin(aj) - C\ loge| ^T[u\. 



The second term in the r.h.s. of (4.79) is simply bounded below as follows 



/ drCoutXini^curl(m,Vu) > - / df eout|i^||Vu|' 
J A J A 



(4.81) 



(4.82) 



We now estimate the third term in the r.h.s. of (4.79): We integrate by parts back to get 



but 



dr Xouti^curl(z«, V«) > - / dr^^ixoutF) ■ (m, Vu) ~ C da \F{1)\ |(m,9,w)| , (4.83) 

A J A JOB 



(4.84) 



drV^iXoutF) ■ {lu, Vu) ^ j dr \^FV ^ {xont) ■ (««, Vm) + 2xout9^B ■ {lu, Vu)} 
and the second term can be bounded using the same computations as in the proof of Lemma 4.2 
2 / dfxontg^B-{iu,Vu) < S [ dfxoutg^\yu\^ + CS-' [ drXoutg^B^\u\^, 

J A J A J A 



where (5 is a parameter to be fixed later. For the first term in (4.84) we use (4.22): 



drV^{Xoni)F-{iu,Vu) 



<Ce-^ /df |V^Xout||r-i?<|g'|u||Vw|< 

J A 



drq^\\Iu? 



{Vxout^O} 



C_ 

fe2 



{Vxout^O} 



drg^\u\ 
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The second inequality uses |1 — i?<| ex e\ logej and (4.69). 

We inject the preceding computations in (4.83), taking into account that B < Ce^^. We also note that 



we have Xout 7^ and/or Vxout 7^ only in the unpleasant set and the average set, so 
/ dfxout-Fcurl(m, Vu) > 



CS 



dfg^\Vu\ 



UPScUASc, 



C_ 

fe2 



UPSa^UASc. 



drg^\u\^-C da \F{l)\\{iu,dru)\ . (4.85) 



We gather equations (4.77), (4.78), (4.79), (4.81), (4.82) and (4.85) to obtain (recalling that |xin| < 1) 



/ dfjg^ |Vm|^ + Fcurl(m,Vii)| 



> 



27r^Cin(aj)|dj| 



(l-7)(i-.).^(a.)|log.|(l-ci^ 



\na,)\ 



(I-7) / dfUt9"\yu\'- / drUt\F\\^u\'-C da\F{l)\\iiu,dru)\ + 
JA JA JdB 

(j-S) f drg^lWuf ^^ f df5>|2 - C\ logep^^M. (4., 

JA "£ JUPS^UASa 



We now choose the parameters in (4.86) as follows: 



7 = 2(5 = 



log I log t 
|loge| 



-log I log e| 
a = a- 



|loge| 



(4.87) 



where a is a large enough constant (see below). This choice allows to bound the terms in (4.86) from 
below: Indeed, if ilo > 2(37r)^^, we have from Proposition A. 2 



i5'(a,)|loge|-|i^(a,)|> 



C 



el logeP 



for any dj G A and thus 



1 



(l-l) --a g'ia,)\loge\ 1-C 



,l0g|l0g£| 

|loge| 



-|^K)I> 



>,,.„„.,, (,-c^)-,n.)i.^>o ,4., 



where we have used (2.25) 



On the other hand, by the definition of ^out, for any re supp(^out), we have either 

|r-i?<|<Ce|log£|-i 

or 

\r~l\ <Ce\loge\'\ 

Therefore, using (4.22), we have in the first case 



(4.89) 
(4.90) 



\F{r)\<C. 



log el 
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In the second case (4.22) yields 



but (2.30) shows that, if f satisfies (4.90), 



\F{r)\ < 



g'ir) > 



C 



e|loE£p 



C 



We conclude that 



e|loge| 
f{r)>C\\oge\\F{r)\:^ \F{r)\ 



for any re supp(^out) and thus 

(1 - 7) / df Coutff'IVtip - / dfCout|F||Vw|2 > 0. 

JA J A 



(4.91) 
(4.92) 



Finally we have from pr86| , (|4^ and pr92| ) 

df (52 |Vu|' + F curl(m, Vu)| > C^^^^^ 
logs 



/ drg2|yy|2_ 



c 



£2 1og|loge| JuPS^UAS. 



f dfg^\u\^-C f da\F{l)\\{iu,dru)\-C\\oge\-'^J'[u]. (4.93) 

JuPScUASa, JdB 



Adding 



A 



dr^{l-\u\^Y- d<jF{l){iu,dru) 



to both sides of (|4.93|) and using (|4.76|), we get the lower bound 

Ifg \u 



s[u]>c\^^:f[u]- Ii°^'" 



|loge| 



£2 1og|loge| JuPS^UAS, 



J^ 2| |2 

dr " '"■' 



da F{l){iu,dru) 



(4.94) 



valid for e small enough and flo > 2(37r) ^. But g^\u\'^ = jvI/'^Pp < Ce ^|log£| ^, whereas the side 
length of a ceU is ©(£1 log£|), thus 



C 

dfg^\u\^ < — j r 

UPS^uAS^ e|log£| 



\UPS^ UASa.\< Cs\ log£| (TVUP + iV^) 



Using (|430|), (|4i36|) and (|4l37|), we deduce 

dr5>p < C£|log£|iV^ < Ce^e^Tiu]. 



(4.95) 



(4.96) 



;7PS„UASo 



On the other hand (2.30) implies that g^{l) > Ce ^\ log£| ^. Combining this fact with the upper bound 
( [2?2l ) yields 

lul <C ondB 



and thus, using (4.21) and Cauchy-Schwarz inequality, 

da F{l){iu,dru) 



dB 



<C{ I da\dru\ 

dB 



1/2 



Using Lemma 4.4 we conclude 

/ daF{l){iu,dri 

JdB 



<c(\loge\'^'\£[u]\' 



/2 



^-N-)^ 



(4.97) 



(4.98) 



(4.99) 
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Combining (4.941, (4.96) and (4.991, we have 



I |log£| log I log e| eV4 



1/2 



(4.100) 



Recall the choice of a in (4.87): We choose now a constant a > 2. Then 

|loge|i^" 



|log£| ^ 

log I log e I log I log e 

and thus there exists a finite constant c such that 

/log I log £ 



< 



log I log e I 
|log£| 



0{e°°) > £[u] > c 



V llog^l 



^M-|iog£|V2|5[„]|i/2_M£^^[,]i/2 



(4.f01) 



where the upper bound comes from (4.4) 



Since the sign of £[u] is not known, we might have two possible cases: If £[u] > 0, (4.4 1 implies that 
\£[u]\ < 0{e°°), which can be plugged in (4.101) yielding 



Oien > £[u] > c ('i^^l^^M 



I log e| 1/4 



|l0g£| 



rl/4 



^MV2 



(4.102) 



This implies 



J'iu] < C— 



|log£|5/2 



el/2 log I ioge|2 



which concludes the proof of Proposition 4.1 if £[u] > 0. 
On the opposite, if £[u] < 0, either 

n^pri_L 11 |i/2|pr i|i/2 ^ f \og\\oge\ |log£|i/4 
0>£:[u]+c|log£| / \£[u]\ ' >c{ T[u] -^j^J'ii 

which implies the result as before, or \£[u]\ < C\ log£|, which gives 



,1/2 



^11 |3/2 ^ /^l0g|l0g£| ^. , |l0g£|l/'' 

' ^ ' - " (, |iog£| -^[""J ^IT^-^M 



1/2 



and thus again (4.8) and (4.9). 



n 



As already noted, the end of the proof could be formulated as an induction. Plugging the estimates 



of Lemma 4.2 in (4.100) and using the upper bound on £[u\ would yield improved estimates of J^[u\ and 



|£[m]|, thus reducing the number of bad cells and improving the boundary estimate. The process could 
then be repeated a large number of times, proving that there are no bad cells at all. The second term in 



the lower bound (4.100 ) would then vanish and the process stop when the first term would reach the order 



of magnitude of the last one (coming from the boundary estimate), thus giving the results of Proposition 

EH 



5 Energy Asymptotics and Absence of Vortices 

In this section we conclude the proofs of our main results. The proof of the energy asymptotics is a 
straightforward combination of the results of Sections [3] and |4j 
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Proof of Theorem 
From Proposition 



which reduces to 



n 



3.1| we have, using the simphfied notation of Section l4l 



pGP > £,GP _ r UogeP^" 

- ^'"« £1/2 log I log £| 



thanks to (4.9 1. Using a regularization of gA.uio ^^ ^ trial function for the functional £^^ as in the proof 



of Proposition 3.1 we get E^^ ^ > E^^ -0{s°^) and thus 

|loge|3/2 



- "" £l/2 1og|log£| 



We conclude the proof of the lower bound recalling that, by definition. 



For the upper bound 



^GP=inf^GP<^^P 



^GP < ^GP ^ ^GP 



GP 



we simply use g,^^ ^(r) exp{i([r2] — Wopt)^} as a trial function for £ 

The proof of the absence of vortices requires an additional ingredient: 

Lemma 5.1 (Estimate for the gradient of u^a)- 

Recall the definition of u^^ in (3.5). There is a finite constant C such that 



n 



iVu. 



'i^o\\L^(A) 



A.<C 



loge 



3/2 



r3/2 



(5.1) 



Proof. We use the short-hand notation defined at the beginning of Section H (in particular u = u^^). 
Recall the variational equation (4.271 



- V (g'^Vu) - 2ig^B ■ Vu + ^/ (|up - l) u = Xg^u. 
From this equation we get the pointwise estimate 

\Au\ < C f ^|V7i| + \B\\\/u\ + \ |g2 (|u|2 -l)u\ + |A| \u\ 
\ g e^ I ^ 

holding true on A. Recalhng that j^I'^^l < e^^/^j logej^i/^ and 

C 



cr-L 



on A, 



we have 



Thus, using (2.25), 



eV2| ioge|3/2 
u\ < C|loge| on J^ 



1 I 2/1 ,2 i\ I ^ ^llogeP 
^\g [\u\ -l)u| <C^^^. 



(5.2) 
(5.3) 

(5.4) 

(5.5) 
(5.6) 
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On the other hand, estimating the chemical potential with ( 4.28 ) and plugging in the results of Proposition 
4.11 we have 



on A. 

Combining p^, (p^, (|4J1) and ^U^ with ([SJ]) and ((5J1), we deduce from Q 

, llogep 



(5.7) 



lAul 



L~(^) 



.-l^ll-l 



L^{A) ^ £3 



From Gagliardo-Nirenberg inequality [N] Theorem at pg. 125], we deduce 



\\Au 



<^, >g^||. ||l/2 II ||l/2 



I log el 



(5.8) 



Inserting (5.5), we conclude 



\Au\ 



L~(^) 



<C 



loge|' 



and we get (5.1) by using (5.5) and the Gagliardo-Nirenberg inequality again. 



n 



We are now in position to complete the 
Proof of Theorem \1.1\ 



The proof relies on a combination of (4.81 and (5.11, as in [BBHlj . 
Suppose that at some point tq such that 



we have 



R^+y\loge\-' <ro<l. 



\\u{f,)\ -l\>e'/'\logef 



Then, using (5.1), there is a constant C such that, for any re Bq := ;B(ro7 C'£^'^/*| logep/^), we have 

||"(rO|-l|>i£^/*^|loge|^ 

This implies (recall ( |4.7[ )) 

,2.2^C|loge|3 



f A-9 f^ I |2^2 ^ C\\o 



/2 



and thus (note that by the initial condition on tq. Bo C yl) 

C|loge|3 



F\u\ > 



rl/2 



(5.9) 



which is a contradiction with (4.8 1 



We have thus proven that (recall (5.5 1 ) 



I* 



GP|2 „2| ^ „2 |L,|2 



g' <9' M-i <C 



log el 

.7/8 



(5.10) 



on A = {-Rh + £| loge| ^ < r < l}. The result then follows by a combination of (2.30) and (5.101. D 
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Remark 5.1 (Absence of vortices in a larger domain). 

By direct inspection of the proof of Theorem |1.1[ one can easily reahze that we could have proven the 

main result, i.e., the absence of vortices, in a domain larger than A, i.e., there is some freedom in the 

choice of the bulk of the condensate. 

More precisely the choice of a larger domain would have implied a worse lower bound on g^ via (2.301 



and in turn a worse remainder in (5.10), but at the same time this would have allowed the extension of 



the no vortex result up to a distance of order e\ loge| " from _Rh for some power a > 1. 
We have however chosen to state the main result in A for the sake of simplicity. 

The proof of the result about the degree of vJ/GP jg g^ corollary of the main result proven above: 
Proof of Theorem |_?.5| 



We first note that the pointwise estimate in (5.10) implies that \1>*^p does not vanish on dB, so that its 
degree is indeed well defined. We then compute 



2.deg{*«^aS} = -^jJa L^9. (^^^ 



do- Or 

u 



dB 

27r([0] -wo)-^ 






-i([n]-uia)^ 



as 



do- Or 

U 



Then 



dB U 



< 



da 



dB 



dr 



<C I dcr \drU\ , 
'dB 



where we have used that \u\ is bounded below by a constant on dB. 

Finally, combining (4.39) and the results of Proposition 4.1 we obtain (recall that g^ > Ce' 

on dB) 



I Aa \dru\ 

JdB 



< 



C|l0g£|^ 

e(log|loge|)2 



Using the Cauchy-Schwarz inequality, we thus conclude from (5.11), (5.12) and (5.13) that 
deg{*^P, dB} ^[n]-ijo + (e^^^^l logep/^^Qg | ^^^^^y 



(5.11) 



(5.12) 



log el 



(5.13) 



(5.14) 



There only remains to recall that (see (3.10)) 

2 

ujo = 



3^/^^e 



e)(e-i|log£|-i/2) 



and that an identical estimate applies to Wopt (see (3.16)) 



Remark 5.2 (Degree of a GP minimizer) . 



U 



According to (5.14), we could have stated the result (1.21) about the degree of "^^^ in terms of 



Wo, I.e. 



the optimal giant vortex phase when the minimization problem is restricted to the annulus A, instead 



of Wopt, i-e., the optimal giant vortex phase in the whole of B. Moreover the remainder in (5.14), i.e., 
0(e~^/^| log£p/^(log I loge|)^^), is much better than the one contained in the final result (1.21), i.e., 
©(e"^! loge|~^/^), which is inherited from (3.10) and (3.16). Note however that the latter remainder 



is the best precision to which one can estimate the giant vortex phase in terms of the explicit quantity 
2/(3-v/7r)e~^. For this reason and the fact that Wopt occurs more naturally in the analysis, we have used 
it in (|1.21|. 
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Appendix A 

In this Appendix we discuss some useful properties of the TF-like functionals involved in the analysis as 
well as the critical angular velocity flc for the emergence of the giant vortex phase. 

A.l The TF Functionals 

We start by considering the TF functional defined in ( |1.9[ ): 

£^^[p]:= /df {-flV^p + e-V}. 
Jb 

Its minimizer over positive functions in L^{B) is unique and is given by the radial density 



1 f-2o2 



pTF(^).^^[^2^TF^^2^2^2] 



[r'-B^l 



(A.l) 



where [ ■ ]+ stands for the positive part and the chemical potential is fixed by normalizing p^^ in L^{B), 
i.e., 

-n^Rl. (A.2) 



,,TF pTF I ^-2 II TF | 

fj, . — £j + £ Up 



Note that, if fJ > 2/(-\/7re), the TF minimizer is a compactly supported function, since it vanishes outside 
A^^ , i.e., for r < i?h, where 



Ry. :— . 1 — 



2 
/neft 



The corresponding ground state energy can be explicitly evaluated and is given by 



i^TF ^-n^ii- 



3^/nefl 



(A.3) 



(A.4) 



By (A.3) and (A.4) the annulus A^^ has a shrinking width of order e| logej and the leading order term in 



the ground state energy asymptotics is —fl^, which is due to the convergence of p'^^ to a delta function 
supported at the boundary of the trap. 

In other sections of the paper we often consider the restrictions of the functionals to domains V strictly 
contained inside B (denoted by f-p^). However in the case of the TF functional there is no need to make 
a distinction between S"^^ and £^^ since all the ground state properties are basically independent of the 
integration domain, provided .4"^^ C T). 



Another important TF-like functional is defined in (1.26) and includes the giant vortex energy con- 
tribution, i.e.. 



Jb 



-2„2 



P'} 



df 



{{[n]-co)' 



r-^p + e-^p^\~2n[n^uj], (A.5) 



where the potential B,^ is defined in (1.25), a; € Z and we have used the normalization in L^{B) of the 
density in the last term. 



The minimization is essentially the same as for (|1.9|): The normalized minimizer is 

P^ir) := 



2„-2 



p^'^-m-ojyr 



and the normalization condition becomes 



logRl 



TTe^{[n]-ujy 



(A.6) 



(A.7) 
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where we have denoted 



p2 im-^) 



M, 



TF 



(A., 



With such a definition the minimizer (A.6) can be rewritten in a form very close to the TF minimizer 



(A.l|,i.e., 



pL' (r) 






2 r 



R' 



In order to mafce a comparison it would then be useful to evaluate the radius R,^ but the equation (A. 7 1 



has no explicit solution. However, since the right hand side of (A.7) tends to zero as e — >■ 0, we can 
expand the left hand side assuming R^^ = 1 + S for some (5 <C 1: 



5^^-5^' + ^('"> = i?(PF^' 



which yields 
1 

We thus have 



-1 = 5 = 



/^e{[n] ~ uj) 



2 1 

3V^e{[n] -u) 97re2([r)] - w)2 ^ ^ [ 



and whether i?^ is larger or smaller than i?h depends in a crucial way on the phase uj: In particular in 



the case of the giant vortex phase ujq (see Proposition 3.2), the sum of the two last terms in the above 
expression vanishes to the leading order (see (3.10)), i.e., it is much smaller than 0{e'^\ logep). 
The ground state energy E"^^ is easy to compute: 



E^' ^~2nm-oj) 



7re2([0]-a;)4 



-fi2 






4a; 



2n{[n]-n) + oie-^\\oge\-^), 



3^/7re 3v^£ 37r£2 

and, assuming that |a;| < 0{e^^), one can easily recognize that the leading term and the first remainder 



coincide with (A.4), i.e., the energy E^ is equal to £' up to second order corrections: 



^TF ^ ^TF 



3V7T£ 



97r£2 



- 2f^ {[n] -n)+ ©(e-^i ioge|-2). 



(A.IO) 



This formula implies that E"^^ is minimized by a phase which is given up to corrections of order 

e-i|loge|-i by 

2 (A.ll) 



c.TF:= 



iy/ne' 



and the same is true for the giant vortex phases Wq (see Proposition 3.2 and (3.10l) and Wopt (see 
Proposition 3.3 and (3.16)). 
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A. 2 The Critical Angular Velocity and the Vortex Energy 

The last part of this Appendix is devoted to the study of the critical velocity ilc, which is defined as the 
angular velocity at which vortices disappear from the bulk of the condensate. To estimate this velocity, 
according to the discussion in SectionWl we have to compare the vortex energy cost ^g'^{r)\ logsj and the 
vortex energy gain |F(r)| (see ( 4.10[ )): This leads to the definition of the function 

H{r) := hXu„{r)noge\ - \F^,{r)\ , (A.12) 

which yields the overall energy contribution of a vortex at a radius r inside the bulk: If H is positive in 
some region, then a vortex is energetically unfavorable there, and, if this holds true in the whole of the 
bulk, the condensate is in the giant vortex phase. 

Before studying the behavior of the above function H, it is however convenient to obtain an explicit 
approximative value for the critical velocity and to this purpose we replace the density g^ with p^^ and 
study the function 

ijTF(r) :=i|loge|pTF(,)_|FTF(,)|, (A.13) 

where the cost function F"^^ is explicitly given by 



FTF(r) := 2 [' ds B^^.{s) ■ e^p^^(s) = e^n^ f ds [^s - [p] - c^^^) s'^] 
JRb JRb 



^TF^,-l](,2_^2)_ (A.14) 



with w'^^ defined in (lA.ll) 



In order to investigate the behavior of the infimum of H"^^ inside the bulk, it is convenient to rescale the 
quantities and set 

z:=en{r^ -Rl), (A.15) 



so that z varies on a scale of order one, i.e., more precisely z G [0, S/^/tt] (see (A. 3)). With such a choice 
the gain function can be easily estimated: 



^TF(^) = '-^1^ """^ dt t [n{t + Rl) M + ^T^] (i + Riy' 



2^2 rr'-Rl 



dtti nt 



3^/7^£ 



o(i; 



^eil 



+ t 



1 

2e 



37^ + ''^'^){'-V^ + ^ 



-— ds s { s — ^ _ 
2e Jo \ 3^ 



0{\\oge\) = 



6e 



where we have used the approximation [1 — 0{{efl) ^)] ^ = 1 + 0{{eft) ^). 
Applying the same rescaling to the energy cost function, we thus obtain 



where 



ir^^(z) =ma-2z 



H'Hr) 



+ 0{e\\oge\) 



zWjz) 
12e '■ 



= 3r2o - 2z 



-z\-0{e\\oge\ 



since z < ^/y/n by the definition of the scaling. Now it is very easy to see that 

H'^^iz) > 3f^o - 2tt-^ - 0{e\ loge|). 
The above considerations lead to the following 



0{\\oge\), 



(A.16) 



(A.17) 



(A.18) 
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Proposition A.l (TF vortex energy). 

For any Q,o > 2(37r)^^ and e small enough, there exists a finite constant C such that 

H^^{r)>Ce-^\\oge\-^ >0 
for any f such that r > _Rh + £| loge|~^. 



Proof. It is sufficient to collect (A. 16), (A. 18) and recall the rescaling (A. 15 1 



n 

We now go back to the original function H and show that the above property holds true as well, i.e., 
r^c = 2(37r)~^e~^| logej"-'^ is the critical velocity for the disappearance of vortices from the bulk of the 
condensate: 

Proposition A. 2 (Critical angular velocity). 

// ilo > 2(37r)^^ and e is small enough, there exists a finite constant C such that 

H{r) >Ce-^\\oge\-^ >0 
for any f such that r > _Rh + £| loge]^^- 
Proof. The result basically follows from what is proven about H^^: We are going to show that 



sup|ijTF(r)-iJ(r)| <Ce-^\logs\-\ 

re A 



In order to prove the above inequality, we use the estimates (2.30) and (2.25) to get 



sup \p^^ir)-gX.M < Ce'/'\loge\y' \\p^^\\^ < Ce-^^l logs|^/^ 
feA 



(A.19) 



(A.20) 



and 



sup|F^^(r)-F„,(r)| <2 / ds \B^ois)\ 9Xu.o(')+ 

re A 



R> 
R< 



Ce il Wo - uj 



TFI 



/ dss-\s'^Rl) + 2snp\p''^{r)~gX^^{r)\ f ds \B, 
jRh re A JRh 



,(5)1 < 



C 



e-' \R^ - R^laX^^^iRy) + e^n'\loge\ + e-'^'\loge\'^' 



< 



C[|logerV^^(^>)+e-^|logeri] <Ce-'\loge\-\ (A21) 



where we have used (A.20), the monotonicity of gA.ujoi''") i^'^^ Propos ition 2.3) and the estimate (3.10). 
Hence one obtains (A.19) and the final result follows from Proposition A.l if r > i?h + e| logej^^^. Indeed, 
using (lA.lSl), (lA.iel) and (lA.lSl) we have 



H^''ir)>Ce-'\loge\-'/^ 

in this case. 

On the other hand, if Ru + s\ logel^^ < r < R^, + e\ logel^/^ it follows from ( |4.22[ ) that 

1 



H{r) > i| log£|5i^„(r)(l - C\ loger^) > C 



e|loge|' 



where the last inequality comes from (4.71 



D 
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